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Abstract 

The full LEP-1 data set collected with the ALEPH detector at the Z pole dur- 
ing 1991-1995 is analysed in order to measure the r decay branching fractions. The 
analysis follows the global method used in the published study based on 1991-1993 
data, but several improvements are introduced, especially concerning the treatment 
of photons and vr'^'s. Extensive systematic studies are performed, in order to match 
the large statistics of the data sample corresponding to over 300 000 measured and 
identified r decays. Branching fractions are obtained for the two leptonic chan- 
nels and eleven hadronic channels defined by their respective numbers of charged 
particles and vr'^'s. Using previously published ALEPH results on final states with 
charged and neutral kaons, corrections are applied to the hadronic channels to derive 
branching ratios for exclusive final states without kaons. Thus the analyses of the 
full LEP-1 ALEPH data are combined to yield a complete description of r decays, 
encompassing 22 non-strange and 11 strange hadronic modes. Some physics impli- 
cations of the results are given, in particular related to universality in the leptonic 
charged weak current, isospin invariance in oi decays, and the separation of vector 
and axial-vector components of the total hadronic rate. Finally, spectral functions 
are determined for the dominant hadronic modes and updates are given for several 
analyses. These include: tests of isospin invariance between the weak charged and 
electromagnetic hadronic currents, fits of the p resonance lineshape, and a QCD 
analysis of the nonstrange hadronic decays using spectral moments, yielding the 
value as{ml) = 0.340 ± O.OOScxp ± 0.014th. The evolution to the Z mass scale yields 
as{M'^) = 0.1209 lb 0.0018. This value agrees well with the direct determination 
from the Z width and provides the most accurate test to date of asymptotic freedom 
in the QCD gauge theory. 
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1 Introduction 



Because of its relatively large mass and the simplicity of its decay mechanism, the r 
lepton offers many interesting, and sometimes unique, possibilities for testing the Standard 
Model. These studies involve the leptonic and hadronic sectors and encompass a wide 
range of topics, from the measurement of the leptonic couplings in the weak charged 
current, providing precise universality tests, to a complete investigation of hadronic 
production from the QCD vacuum. In the latter case, the r decay results have proven to 
be complementary to those from e^e~ data, enabling detailed studies at the fundamental 
level to be performed, through the determination of the spectral functions which embody 
both the rich hadronic structure seen at low energy and the quark behaviour relevant 
in the higher energy regime. The spectral functions play an important role in the 
understanding of hadronic dynamics in the intermediate energy range and constitute 
the basic input for QCD studies and for evaluating contributions from hadronic vacuum 
polarization. The latter are needed for precision tests of the electroweak theory through 
the running of a to the Mz scale and to compute the muon anomalous magnetic moment. 
The topics of interest include: testing the isospin invariance between the weak charged 
and electromagnetic hadronic currents; evaluation of chiral sum rules making use of the 
separate determination of vector and axial- vector components; and a global QCD analysis 
including perturbative and nonperturbative contributions, offering the possibility of a 
precise extraction of the strong coupling at a relatively low energy scale, thus providing 
a sensitive test of the running when compared to the value obtained at the Mz scale. 

A global analysis of the r branching ratios was performed by the ALEPH Collaboration 
using the data recorded at LEP through the process e+e^ ^ Z ^ r^r^ in 1990-1991 [1] 
and in 1991-1993 [2, 3]. Specific studies of final states including kaons were published 
using an early data set [4, 5], then with the full 1991-1995 statistics [6, 7, 8], and finally 
compiled [9]. Nonstrange spectral functions were determined earlier [10], then separately 
for the vector [11] and the axial-vector [12] components from the 1991-1994 data, leading 
to a precise determination of asirri^) and showing that nonperturbative terms are small 
even at this relatively low scale. The strange spectral function was likewise extracted, 
from which it was possible to determine the strange quark mass [9]. Finally, a specific 
study of T) and cu resonance production in r decays was carried out [13] . 

In this report a complete and final analysis of r decays is presented using the global 
method. All data recorded at LEP-I from 1991 to 1995 with the ALEPH detector are 

used, thus providing an update of the previous results which were based on partial 
data sets. The increase in statistics — the full sample corresponds to about 2.5 times 
the luminosity used in the last published global analysis [2, 3] — not only allows for 
a reduction of the dominant statistical errors but, more importantly, provides a way 
to better study possible systematic effects and to eventually correct for them. Several 
improvements of the method have been introduced in order to achieve a better control 
over the most relevant systematic uncertainties: simulation-independent measurement 
of the rr selection efficiency; improved photon identification especially at low energy, 
where the separation between photons from r decays and fake photons from fluctuations 
in hadronic or electromagnetic showers is delicate; a new method to correct the Monte 
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Carlo simulation for the rate of fake photons; and stricter criteria for channels with 
small branching fractions. For consistency and in order to maximally profit from the 
improved analysis all data sets recorded from 1991 to 1995 have been reprocessed with 
the latest version of the reconstruction and selection programs. The results presented 
in this report thus supersede those already published in Ref. [1, 2, 3]. In the same way 
the new determination of the nonstrange spectral functions replaces earlier results in 
Ref. [10, 11, 12]. Only the measurements on final states containing kaons, which were 
already based on the full statistics, remain unchanged [6, 7, 8, 9]. 

The new analysis has been carried out independently for the two data sets taken in 
1991-1993 and 1994-1995, respectively. Each data set is associated to its own Monte Carlo 
sample describing the corresponding detector performance. This method was chosen as it 
provides an easy way to compare with the previous analysis of the 1991-1993 data sample 
and to check the consistency between the two sets of new results. The latter comparison is 
particularly meaningful as the systematic uncertainties are estimated from the agreement 
between data and Monte Carlo samples, once corrections are applied, and are thus mostly 
uncorrelated. The final results are obtained from the combination of both sets. 

The paper is organized as follows. After a brief description of the ALEPH detector, 
and of the data and simulated samples in Section 2, the main analysis tools are presented 
in Sections 3 through 9: selection of rr events and determination of non-r backgrounds, 
charged particle identification, separation of genuine and fake photons, 7r° reconstruction, 
decay classification, channel-by-channel adjustment of the number of fake photons in the 
simulation, and branching ratio calculation. The determination of the different sources 
of systematic uncertainties is the subject of Section 10, while Section 11 provides some 
global checks on the results. Corrections are applied to raw results in order to obtain 
branching ratios for exclusive channels, as discussed in Section 12. The consistency 
between the results obtained from the two data sets is examined in Section 13 and the final 
combined results are given. Section 14 deals with a few physics topics directly related to 
the measurement of the branching ratios, such as imiversality in the leptonic charged weak 
current, the branching ratio into in the context of vacuum polarization calculations, 
ai decay fractions into the two isospin-related modes Stt and 7r27r°, and separation of 
vector and axial-vector modes. The subject of Section 15 is the measurement of hadronic 
spectral functions, the vector part of which is compared to corresponding results from 
e^e~ annihilation in Section 16, providing tests of the Conserved Vector Current (CVC) 
hypothesis. Finally, a QCD analysis of the vector and axial-vector spectral functions is 
performed in Sections 17 with a determination of the strong coupling at the r mass scale, 
as{m1) and of the small nonperturbative contributions. 

2 The ALEPH detector and the data sample 

A detailed description of the ALEPH detector and its performance can be found 
elsewhere [14, 15]. Only the features relevant to this analysis are briefiy mentioned here. 

Charged particles are measured by means of three detectors. The closest detector to 
the interaction point is a silicon vertex detector ( VDET) , which consists of two concentric 
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barrels of double-sided microstrip silicon detectors. An inner tracking chamber (ITC), 
with eight drift chamber layers, surrounds the vertex detector. The ITC is followed by a 
time projection chamber (TPC), a cylindrical three-dimensional imaging drift chamber, 
providing up to 21 space points for charged particles, and up to 338 measurements of 
the ionization loss, dE/dx. Combining the coordinate measurements of these detectors, 
a momentum resolution Spx/Pr — ^ ' ® ^ " 10~^/pr (with pt in GeV) is achieved 
in the presence of a 1.5 Tesla magnetic field. 

The electromagnetic calorimeter (ECAL), located inside the coil, is constructed from 
45 layers of lead interleaved with proportional wire chambers. The position and energy 
of electromagnetic showers are measured using cathode pads subtending a solid angle 
of 0.9° X 0.9° and connected internally to form projective towers. Each tower is read 
out in three segments with a depth of 4, 9 and 9 radiation lengths, yielding an energy 
resolution SE/E = 18%/\/E + 0.9% (with E in GeV) for isolated photons and electrons. 
The inactive zones of this detector, referred to as cracks, represent 2 % in the barrel 
and 6 % in the endcap regions. The analysis of the hadronic r decays presented in this 
paper benefits from the fine granularity and from the longitudinal segmentation of the 
calorimeter, which play a crucial role in the photon and 7r° reconstruction, and in the 
rejection of fake photons. 

The hadron calorimeter (HCAL) is composed of the iron of the magnet return yoke 
interleaved with 23 layers of streamer tubes giving a digital hit pattern and has a projective 
tower cathode pad readout of hadronic energy with a resolution of about 85%/ V^- 
Outside this calorimeter structure are located two additional double layers of streamer 
tubes, providing three-dimensional coordinates for particles passing through the HCAL. 

The trigger efficiency for r pair events is measured by comparing redundant and 
independent triggers involving the tracking detectors and the calorimeters. The measured 
trigger efficiency is better than 99.99% within the selection cuts. 

Tau-pair events arc simulated by means of a Monte Carlo program, KORALZ07, which 
includes initial state radiation computed up to order and exponentiated, and final 
state radiative corrections to order a [17]. The simulation of the subsequent r decays also 
includes single photon radiation for the decays with up to three hadrons in the final state. 
Longitudinal spin correlations are taken into account [18] and the sin^ 9^ value is adjusted 
to be consistent with the determination from the measurement of r polarization [16] . The 
GEANT-bascd simulation, with the detector acceptance and resolution effects, is used to 
initially evaluate the corresponding relative efficiencies and backgrounds. It also includes 
the tracking, the secondary interactions of hadrons, bremsstrahlung and conversions. 
Electromagnetic showers are generated in ECAL according to parameterizations obtained 
from test beam data [14]. For all these effects, detailed comparisons with relevant data 
distributions are performed and corrections to the MC-determined efficiencies are derived, 
as discussed below for each specific problem. 

The data used in this analysis have been recorded at LEP-1 in 1991-1995. The numbers 
of detected r decays are correspondingly 132 316 in 1991-1993 and 194 832 in 1994-1995, 
for a total of about 3.3 • 10^. The ratios between Monte Carlo and data statistics are 7.3 
and 9.7 for the 1991-1993 and 1994-1995 periods, respectively. Monte Carlo samples were 
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generated for each year of data taking in order to follow as closely as possible the status 
of the detector components. For convenience and easier reference to previously published 
results, the total data and Monte Carlo sets corresponding to the 1991-1993 and 1994- 
1995 periods are considered separately. Also, part of the data (approximately 11% of the 
r sample) was taken off the Z peak in order to measure the lineshape. Corresponding 
Monte Carlo sets were generated to account for these conditions. 

3 Selection of rr events 
3.1 The rr event selector 

The principal characteristics of rr events in e^e~ annihilation at the Z energy are 
low multiplicity, back-to-back topology and missing energy. Particles in each event are 
reconstructed with an energy flow algorithm [15] which calculates the visible energy, 
avoiding double-counting between the TPC and the calorimeter information. The 
reconstructed event is divided into two hemispheres by a plane perpendicular to the 
thrust axis. The jet in a given hemisphere is defined by summing the four-momenta of 
energy flow objects (charged and neutral). The energies in the two hemispheres including 
the energies of photons from final state radiation. Ex and £^2, are useful variables for 
separating Bhabha, /i/i and 77-induced events from the rr sample, while relatively larger 
jet masses, wider opening angles, and higher multiplicities indicate Z ^ qq events. 

All these features are incorporated in a standard selector used extensively in ALEPH. 
While more detailed information can be found in earlier publications [2, 3, 19] the most 
important cuts are listed below: 

1. The total charged multiplicity should be at least two and no more than eight. 

2. Each hemisphere is required to have at least one charged track. 

3. The scattering angle 9* in the rr rest frame, calculated using the measured polar 
angles 9i and 82 of the two jets through the relation cos^* = sin( ^^~^^ )/ sin(^i^), 
should satisfy | cos^*| < 0.9. 

4. The acoUinearity angle between the two jets should exceed 160°. 

5. The sum of the jet energies should be larger than 0.35 x Ebeam and the difference 
between the transverse momenta of the two jets should satisfy \ Apt\ > 0.066 x Ebeam- 

6. At least one track should extrapolate to the interaction point within ±1 cm 
transversally to and ±5 cm along the beam axis. 

7. A r-like hemisphere is defined by its charged multiplicity equal to one and its 
invariant mass less than 0.8 GeV. For non r-like hemispheres, the product of the 
number of energy fiow objects in each hemisphere should be less than 40 and the 
sum of the maximal opening angle between any two tracks in each hemisphere less 
than 0.25 rad. 
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8. The ratio of the total energy to the beam energy should be smaller than 1.6, except 
for ////-like events (both leading tracks are muons, or one is a muon and the energy in 
the opposite hemisphere is greater than 0.9 x -Bbeam) where this cut is relaxed to 1.8. 
For Bhabha-like events (all the charged tracks are electrons), the cut is tightened 
to 1.4 if the tangent at the origin to the leading electron points to within ±6 cm of 
an ECAL crack. 

9. As an additional cut against the residual Bhabha background, the total energy 
measured independently in the wire planes of ECAL has to be less than 1.8 x Ebeam- 

10. Cuts using tight matching in space and momentum between the tracks in opposite 
hemispheres are used to reject cosmic rays [2]. 

In the previous analysis [2, 3] additional cuts had been introduced in order to further 
reduce the contamination from Bhabha and e^e~ — > /i "*"/*" processes. In the present 
work it was chosen to simplify the procedure in order to conveniently measure selection 
efficiencies on the data, at the expense of a slightly larger background contamination 
which is anyway also measured in the data sample as explained below. 

3.2 Measurement of selection efficiencies 

The "break-mix" method introduced for the determination of the rr cross section [19] is 
used to measure the selection efficiency of the cuts listed above. A selection of rr events is 
implemented using very tight cuts on only one hemisphere. Then the opposite hemisphere 
corresponds to an essentially unbiased r decay, which is stored away. Pairs of selected 
hemispheres are combined to construct a rr event sample built completely from data. 
This sample is used to measure the efficiency of the cuts based only on energy and/or 
topology. 

The same procedure is applied to the Monte Carlo sample as well as to the data and 
the ratio of the corresponding efficiencies is taken as a correction to the pure Monte Carlo 
efficiency. In this way the small correlations between hemispheres which are neglected in 
the break- mix data sample are properly taken into account. 

This method is used to determine the efficiencies of cuts 5 to 9 listed above for both 
MC and data for large branching ratio channels. Channels with small branching ratios are 
corrected using the same factor as for nearby topologies, with an error properly enlarged. 
Results are given in Section 10.2. The measured efficiencies are found to be very close to 
those obtained by the simulation, deviations being at most at the few per mille level. This 
situation stems from the facts that the r decay dynamics is — apart from small branching 
ratio channels — very well known, the selection efficiencies are large and the simulation 
of the detector is adequate. 

The overall selection efficiency of rr events is 78.9%. This value increases to 91.7% 
when the rr angular distribution is restricted to the detector polar acceptance (cut 
3), giving a better indication for the efficiency of the cuts designed to exclude non- 
rr backgrounds. In addition, when expressed relatively to each r decay, the selection 
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efficiencies are only weakly dependent on the final state, with a maximum difference of 
10% amongst the considered decay topologies. 

3.3 Estimation of non-rr backgrounds 

A method — developed for the measurement of the r polarization [16] — has been used 
to measure the contributions from the major non-r backgrounds: Bhabhas, fi'^fi~ pairs, 
and 7*7* —>■ e'^e~ , , and hadrons events. The basic idea is to apply cuts to the 
selected rr data in order to reduce as much as possible the rr population while keeping 
a high efficiency for the background source under study. The procedure does not require 
an absolute normalization from the Monte Carlo simulation of these channels, only a 
qualitative description of the distribution of the discriminating variables. In the following, 
'e' and '/x' designate charged particles identified (see Section 4.2) as electrons and muons, 
respectively. 

For Bhabha background, different event topologies are considered in turn, depending 
on whether one of the charged particle traverses an inactive area of the ECAL between 
detector modules. For e — e and e — crack topologies, either the acoplanarity angle is 
required to be larger than 179° or the acollincarity angle should be more than 175°, in 
which case the difference of transverse energies is required to be less than 3 GeV. The 
corresponding Bhabha efficiency (for the final sample) is 88% as determined from the 
simulation using the UNIBAB generator [20] . The angular distribution of the restricted 
sample is then fitted to rr and Bhabha components (also including a small contribution 
from the other non-r backgrounds) from the simulation. Therefore, the Bhabha Monte 
Carlo input is only used to determine the (large) selection efficiency and the cos 9* shape 
inside the final sample, not relying on any determination of the absolute Monte Carlo 
normalization. The derived Bhabha contribution has a statistical uncertainty which is 
assigned as a systematic error. Several combinations of variables have been tried, showing 
a good stability of the result within its error. Figure 1 illustrates the determination for 
the e — e topology in the 1994-1995 sample. For the more numerous e — hadron sample, 
the same cuts are used, but they have to be supplemented by an additional requirement to 
suppress true hadrons as compared to electrons misidentified as hadrons. This is achieved 
by restricting opposite hadrons to have an electron identification probability larger than 
0.01; most of the true hadrons are below this value. 

A similar technique is used to estimate the yU-pair background. Figure 2 shows the 
corresponding plots for the fi — fi topology, requiring an acoplanarity angle larger than 
178°, for which 90% efficiency for ee — > /x/x events in the final sample is achieved. Here, 
the fitted distribution is that of the calculated photon energy along the beam for a 
postulated ee — > /x//776eam kinematics to take the most general case compatible with 
the only information on the two muons. The ee — > n/j, background signal is clearly seen 
for small energies. 

The remaining background from 7*7* ^fi is easily found in the — /i topology after 
requiring the acoUinearity angle to be smaller than 170° leading to an efficiency of 57%. 
The distribution of the acoplanarity angle shows a distinct signal for the background, as 
seen in Fig. 3. 
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Figure 1: Determination of the remaining Bhabha contribution in the final e-e sample 
for 1994-1995 data: cos 6'* distributions for (a) data, (b) rr simulation, and (c) Bhabha 
simulation. The solid line represents a fit of the data with the relative normalization of 
the TT and Bhabha contributions left free. 
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Figure 2: Determination of the remaining e^e~ — > jJi^ jjT contribution in the final ji-jji 
sample for 1994-1995 data using the kinematically calculated energy of a hypothetical 
photon emitted along one of the momenta: (a) data, (b) rr simulation, and (c) 
e+e" 1^'^ simulation. The solid line represents a fit of the data with the relative 
normalization of the rr and //-pair contributions left free. 
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Figure 3: Determination of the remaining 7*7* contribution in the final 

sample for 1994-1995 data: acoplanarity angle distributions for (a) data, (b) rr simulation, 
and (c) e'^e~ — > {e^e~)^'^ simulation. The solid line represents a fit of the data with 
the relative normalization of the rr and 77-induced //-pair contributions left free. 
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Table 1: Numbers of non-r background events in observed classes as defined in Section 7. 



class 


1991-93 


1994-95 


e 


598 ± 46 


745 ± 58 




409 ± 45 


380 ± 40 


h 


93 ± 11 


100 ± 13 




141 ± 22 


178 ± 26 


/i27r° 


44 ± 9 


81 ± 16 


/i37r° 


26 ± 7 


35 ± 9 




12 ± 3 


19 ± 5 


3h 


87 ± 20 


129 ± 30 




97 ± 23 


165 ± 39 


3/t27r° 


27 ± 7 


36 ± 10 


3/i37r° 


13 ± 4 


25 ± 7 


5/i 


3 ± 1 


7 ± 2 




16 ± 5 


21 ± 6 


Class 14 


249 ± 38 


303 ± 52 


sum 


1815 ± 86 


2224 ± 107 



The same procedure is applied for 7*7* — > ee, but additional cuts (for example 
requiring the total event energy to be less than 35 GeV) must be applied to decrease 
the Bhabha contribution. Similar techniques are used to determine the contribution of 
7*7* — > hadrons for low multiplicity events. 

Cosmic ray background is determined to be negligible [2]. All other contributions 
are estimated with the proper Monte Carlo generators. This is the case for four- 
fermion processes and hadronic Z decays whose contribution is estimated from the Lund 
generator [21]. Tests were previously made to ascertain the reliability of the prediction 
for the rate of surviving low multiplicity events by comparing to data [3], resulting in a 
systematic uncertainty included in the analysis. 

The non-T backgrounds in each channel are listed in Table 1. The contaminations from 
the different sources are given in Table 2 and amount to a total fraction of (1.23 ±0.04)% 
in the full data sample. 

4 Charged particle identification 
4.1 Track definition 

A good track is defined to have a momentum greater than or equal to 0.10 GeV (and not 
smaller than the momentum resolution), | cos(^)| < 0.95, at least 4 hits in the TPC, and 
its minimum distance to the interaction point within 2 cm transversally and 10 cm along 
the beams. If a charged track is not a good one, has at least three hits in the TPC and 
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Table 2: Estimated contaminations from the different background processes. 



process 



contamination (xlO ^) 



Bhabha 
ee — > /i/i 



4.0 
1.0 
2.2 
1.1 
3.9 
0.1 



two-7* processes 

four-fermion 

hadrons 



cosmic rays 



sum 



12.3 



the minimum distance to the interaction point is within 20 cm transversally to and 40 cm 
along the beams, then it is called a bad track. 

In classifying r decays, only good tracks are used, after removing those identified as 
electrons which are used to reconstruct converted photons; the electrons identified as bad 
tracks are also included in reconstruction of conversions. 

4.2 The likelihood identification method 

Charged particle identification plays a crucial role in the measurement of r branching 
ratios. In this analysis, as in the previous ones, a likelihood method is used to incorporate 
the information from the relevant subdetectors. In this way, each charged particle is 
assigned a set of probabilities from which a particle type is chosen. 

Eight discriminating variables are used in the identification procedure: dE/dx in the 
TPC, two estimators (transverse and longitudinal) of the shower profile in ECAL, the 
average shower width measured with the HCAL tubes in the fired planes, the number of 
fired planes among the last ten, the energy measured with HCAL pads, the number of 
hits in the muon chambers, in a road ±4(T-wide around the track extrapolation, where 
(J is the standard deviation expected from multiple scattering, and finally, the average 
distance (in units of the multiple-scattering standard deviation) of the hits from their 
expected position in the muon chambers. 

Probability densities //(a;,) of discriminating variable arc determined using the 
ALEPH simulation for each particle type j, where j = e,fi,h {h = hadron). No attempt 
is made in this analysis to separate kaons from pions in the hadron sample since final 
states containing kaons have been previously studied [6, 7, 8]. Each charged particle is 
assigned to the type with the largest global estimator defined as 



The performance of the particle identification has been studied in detail using control 
samples of Bhabha events, /x/x pairs, 7*7*-induced lepton pairs and hadrons from 7r°- 
tagged T decays over the full angular and momentum range. 
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Figure 4: Hadron identification efficiency and misidentification probability obtained from 
the TT simulation, corrected from data using the control samples, for the 1994-1995 data 
set. 



The measurement of the particle identification efficiencies and the misidentification 
probabilities are described in details in Ref. [2]. Basically the efficiencies are obtained 
from the r Monte Carlo in order to include the relevant detector description, but corrected 
as a function of momentum by the ratio of data over Monte Carlo efficiencies as measured 
with the control samples. This procedure takes care of the slightly different environment 
for a given particle between the various samples. Figures 4 and 5 show the data-corrected 
particle identification efficiencies and misidentification probabilities for e, // and hadron 
tracks as a function of track momentum. Significant differences are observed between data 
and simulation, as in previous analyses, emphasizing the necessity of measuring directly 
these quantities on the data. 



5 Photon identification 

5.1 The problem of fake photons in ECAL 

The high collimation of r decays at LEP energies quite often makes photon reconstruction 
difficult, since these photons are close to one another or close to the showers generated 
by charged hadrons. Of particular relevance is the rejection of fake photons which may 
occur because of hadronic interactions, fiuctuations of electromagnetic showers, or the 
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Figure 5: Lepton identification efficiencies obtained from the rr simulation, corrected 
from data using the control samples, for the 1994-1995 data set. 



overlapping of several showers. These problems reach a tolerable level thanks to the fine 
granularity of ECAL, in both transverse and longitudinal directions, but they nevertheless 
require the development of proper and reliable methods in order to correctly identify 
photon candidates. 

A feeling of the fake photon problem can be obtained from Fig. 6 showing the 
comparison of simulated (at the generator level, after smearing for resolution effects) 
and observed photon energy spectra. A clear difference at low energy indicates a 
large contribution from fake photons. Since the fake photons have a complex origin 
and probably quite dependent of the details of the detector simulation, the procedure 
developed to separate them from the genuine primary photons from r decays needs to be 
calibrated using the data. 



5.2 Likelihood method for calorimeter photons 
5.2.1 Photon reconstruction 

The clustering algorithm for the photon reconstruction [15] starts with a search for local 
maxima among the towers in the three ECAL stacks. The segments of a projective tower 
which share a face in common with the local maximum are linked together into a cluster. 
At the end of the procedure, every segment of a tower is clustered with its neighbour of 
maximal energy. Then a cluster is accepted as a photon candidate if its energy exceeds 
350 MeV and if its barycentre is at least 2 cm from the closest track extrapolation; this 
distance is slightly modified later after corrections for the finite size of the pads are applied. 

The energy of the photon is calculated from the energy deposition in the four central 
towers only when the energy distribution is consistent with the expectation of a single 
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Figure 6: Energy spectra of reconstructed photons in simulated rr events (points) 
and of primary produced photons by KORALZ07 after smearing for energy resolution 
(histogram). The fake photon contribution at lower energies is clearly indicated by 
normalizing to the number of photons with > 5 GeV). The distribution of primary 
photons is obtained from a statistically reduced sample. 

photon, otherwise the sum of the tower energies is taken. 
5.2.2 Method for photon identification 

As presented above, the main problem is the separation between genuine photons from r 
decays and the numerous fake photons generated either by the hadronic tracks interacting 
in ECAL, or by the fluctuation of the electromagnetic showers. The situation is worse for 
low energy photons which are however needed to retain a high 7r° reconstruction efficiency. 

A likelihood method is used for discriminating between genuine and fake photons. For 
every cluster, a photon "probability" is defined 

pgenuine 

p = (2) 

T pgenuine _j_ pfake ' ^ ' 

where is the estimator under the photon hypothesis of type j. By definition, genuine 
photons have a P^ value near one while fake photons have P^ near zero. Each estimator 
Pj is constructed according to 

P' = I[T'i(^i) (3) 

i 

where Vj is the probability density for the photon hypothesis of type j associated to the 
discriminating variable Xi. 
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Figure 7: Genuine-to-fake photon ratio as a function of photon energy from simulation, 
for one- and three-prong r decays, and in barrel and endcap parts of ECAL. 

To take into account the ratio of fake to real photons which is strongly energy- 
dependent, another quantity is defined as 



pgenmne 

P' = - 

7 pgenuine _j_ pfake j p ^ ' 



(4) 



where Rgf is the ratio of the numbers of genuine and fake photons, which depends on the 
photon energy, the number of charged tracks and the number of photons in the hemisphere. 
This quantity is used in 7r° reconstruction in order to reduce the fraction of 7r° candidates 
where one of the photons is fake. Figure 7 shows the genuine-to-fake photon ratio as a 
function of energy as obtained from simulation. 

Discriminating variables for each photon candidate used are the distance to the charged 
track ((ic), the distance to the nearest photon (rf^), a parameter from the clustering process 
resulting from the comparison between the energies in the cluster and in the neighbouring 
cells(G), the fractions of energy deposition in ECAL stacks 1 (-Ri), 2 {R2) and 3 (-R3), and 
a parameter related to the transverse size of the energy distribution (T). The reference 
distributions for photons with energy between 1.0 and 1.5 GeV are shown in Fig. 8. 



5.2.3 Improvements with respect to the previous analysis 

Major improvements were introduced at this stage in the analysis compared to the 
previous one [3]: 
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Figure 8: Reference distributions for fake photons (histogram) and genuine photons 
(points) with energy between 1.0 and 1.5 GeV. The underflows and overflows in each 
plot are also taken into account in the reference distributions. 
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Figure 9: Correlations between Ri and i?2, and between Ri and R3, for photons with 
energy below 3 GeV. 



(1) i?3 is used instead of R2. This change was made because of the large 
correlation between R2 and Ri for genuine photons, whereas the likelihood method ignores 
correlations between the discriminating variables. Also, R3 contributes some additional 
discriminating power against fake photons (see Fig. 9). As a consequence, genuine and 
fake photons are better separated. 

(2) Instead of using the same reference distributions for the whole energy range 
as before, energy-dependent reference distributions are now used to improve the 
discrimination between genuine and fake photons. This is made necessary because the 
reference distributions show pronounced energy dependence at energies below a few GeV's 
(see Fig. 10). 

(3) Large differences are observed between the reference distributions of ECAL barrel 
and endcap photons (Figs. 11 and 12), so different sets of reference distributions are used 
for photons in different regions of the detector. 

By comparing the distributions of photons at high energy, good agreement is observed 
between data and MC photons in dc, d^, and G distributions, while for R^ and T there 
are discrepancies (Fig. 13). The Monte-Carlo-determined reference distributions are then 
corrected using data. In order to do so, the Ri and i?3 data distributions are obtained 
separately for real and fake photons by fitting the uncorrelated dc distributions with fake 
and genuine photon dc distributions from the simulation. For T, the value is shifted and 
smeared according to the difference between data and Monte Carlo. 



21 



0.15 r 




-100 -50 50 100 5 10 15 20 

(cm) (degree) 




Figure 10: Comparison of the reference distributions for genuine photons in the 0.5- 
1.0 GeV (points) and 1.5-2.0 GeV (histogram) energy ranges. 
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Figure 11: Comparison of the reference distributions in the barrel (histogram) and endcap 
(points) parts of ECAL for genuine photons between 1.0 and 1.5 GeV. 
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Figure 12: Comparison of the reference distributions in the barrel (histogram) and endcap 
(points) parts of ECAL for fake photons between 1.0 and 1.5 GeV. 
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Figure 13: Comparison of R3 and T distributions in data (points) and simulation 
(histogram) for photons with energy between 3.0 and 6.0 GeV, in ECAL barrel and 
endcap parts. 

5.2.4 Total set of reference distributions 



Fine granularity is introduced when needed in the reference distributions. Due to the 
observed strong energy dependence of some distributions, more energy bins are considered: 
for genuine photons 7 bins (0.3-0.5, 0.5-1.0, 1.0-1.5, 1.5-2.0, 2.0-2.5, 2.5-4.0, above 4 GeV) 
and for fake photons 4 bins (0.3-0.5, 0.5-1.0, 1.0-2.0, above 2 GeV). Different distributions 
are set up according to the number of charged tracks (1 and > 1) and of photon candidates 
(1, 2, 3, 4, > 4 for one prong and 1, 2, > 2 otherwise). Also the barrel and endcap regions 
of the calorimeter are separated so that the different amount of material in front and the 
different geometry is taken into account. Thus the photon identification procedure relies 
on the construction of 768 reference distributions for fake photons and 1344 for genuine 
photons. 

Figure 14 shows the comparison of the photon identification probability (P^) before 
and after introducing energy-dependent reference distributions and the variable. A 
spectacular improvement is observed in genuine and fake photon discrimination: at low 
energy, a clear contribution of genuine photons can be seen, while at high energy, a small, 
but well identified, fake photon component shows up. Further checking, along the lines 
discussed below in Section 8, shows the validity of the latter signal on the simulated 
sample. 



5.2.5 Photon energy calibration 

Better photon energy calibration is achieved compared to the previous analysis. The 
procedure, aiming at a relative calibration between data and simulation, is implemented 
in several steps. First, a calibration is done using electrons from r decays, treating the 
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Figure 14: Comparison of the data probability distributions of a cluster to be identified as 
a genuine photon before (histogram) and after (error bars) introducing energy- dependent 
reference distributions and the i?3 variable. 
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Figure 15: Reconstructed 7r° mass (upper plot) and resolution (lower plot) after photon 
energy and resolution calibration. Good agreement between data and simulation is 
observed. 



ECAL information as for a neutral cluster. This method is reliable only above 2 GeV 
because the electron track curvature in the magnetic field causes the shower to be ill- 
pointing to the interaction region. The results are then extrapolated to lower energies. 
The bias introduced by the extrapolation is studied with the simulation comparing the 
true and the extrapolated energies, and the same correction is applied to data. Then the 
final calibration is achieved by comparing the reconstructed 7r° mass and resolution as a 
function of energy. The agreement between data and simulation is satisfactory, as seen in 
Fig. 15. 



5.3 Converted photons 

In order to identify photons which convert inside the tracking volume all oppositely 
charged track pairs of a given hemisphere identified as electrons are considered. These 
candidates are required to have an invariant mass smaller than 30 MeV and the minimal 
distance between the two helices in the transverse plane must be smaller than 0.5 cm. 
Finally, all remaining unpaired charged tracks identified as electrons are kept as single 
track photon conversions. These include Compton scatters or asymmetric conversions 
where the other track was either lost or poorly reconstructed. 

The P-y estimator is naturally set to unity for a converted photon. As the conversion 
rate in simulation is found to be lower than that of data by (7.6 ± 1.2) % in 1991-1993 
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Figure 16: Comparison of the distributions of the conversion point radius and the two- 
track invariant mass for photon conversions in data (points) and simulation (histogram) for 
the 1994-1995 data sample. The distribution of material in the detector is well simulated 
except for a global lower rate in the simulation which is already corrected here. 



and (9.3 ± 1.0) % in 1994-1995, some converted photons in data are turned randomly 
into calorimeter photons and assigned a P-^ value, generated according to the observed 
distribution for calorimeter photons at the same energy. Figure 16 shows the comparison 
of the radial conversion point and the invariant mass of two-track conversions between 
data and simulation. After correcting for the conversion rate, good agreement is observed. 



6 TT reconstruction 
6.1 Different tt^ types 

The goal of the tt^ reconstruction procedure is to achieve a high efficiency while keeping the 
"fake" TT^'s at a reasonably small level. Three different kinds of vr'^'s are thus reconstructed: 
resolved 7r° from two-photon pairing, unresolved vr*^ from merged clusters, and residual 
7r° from the remaining single photons after removing radiative, bremsstrahlung and fake 
photons. 
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6.2 Resolved tt s 

Since fake and genuine photons are better separated than previously pubhshed [3], a 
cut on before 7r° reconstruction is introduced to remove a large fraction of fake 
photons, while keeping the real photon efficiency high. The cut is energy-dependent, 
InP^cut > 0.15 - 2.9 X E^, with an efficiency > 90% at < 0.5 GeV, and ~ 100% above 
2 GeV. 

A 7r° estimator Dij is defined [3] to take into account the genuine and fake photon 
ratio at different energies: 

A,- = p; • p; • , (5) 

where Ptt" is the probability of a kinematic 7r°-mass constrained fit of the two photons. 
It should be noted that P^, is used here instead of P^, as in the earlier analysis, in order 
to keep the same confidence for tt^'s reconstructed from photons at different energies. 
Finally, tt*^ candidates are retained if their Dij value is larger than 1 x 10~^. 

In addition, a criterion must be established for choosing among all the accepted i — j 
pairs in a multiphoton environment. Configurations with the maximum number of 7r°'s 
arc retained and among those, the configuration which maximizes the product of all Aj 's 
is kept. 

In the old analysis [3], a significant difference in observed 7r° mass and resolution 
in data and simulation was observed. After the recalibration of the photon energy and 
resolution using the procedure described above, data and simulation agree reasonably well 
as shown in Fig. 15. 

6.3 Unresolved 7r°'s 

As the vr" energy increases it becomes more difficult to resolve the two photons and the 
clustering algorithm may yield a single cluster. The two-dimensional energy distribution in 
the plane perpendicular to the shower direction is examined and energy-weighted moments 
are computed. Assuming only two photons are present, the second moment provides a 
measure of the 77 invariant mass [15]. 

Figure 17 shows this invariant mass distribution for data and simulation. Clusters with 
mass greater than 0.1 GeV are taken as unresolved vr'^'s. Some discrepancy is observed 
between data and Monte Carlo and in order to keep the same efficiency a slightly different 
value for the tt" mass cut is applied in the simulation. It should be noted that this cut 
only affects the definition of unresolved and residual 7r°'s, and thus does not change in 
first order the 7r° multiplicity, hence the definition of the r decay final state. Thus it 
has essentially no influence on the branching ratio analysis, since both unresolved and 
residual 7r°s are used in the r decay classification. 

6.4 Residual tt^'s 

After the pairing and the cluster moment analysis, all the remaining photons inside a 
cone of 30° around the jet axis are considered. They come from different sources: 



29 



> 

0) 



LU 



> 



LU 



> 
CD 



20000 
15000 
10000 
5000 


1500 
1000 
500 

150§ 

1000 

500 



800 

600 
400 
200 




E < 8 GeV 




<E <15GeV 



1 5 < E < 25 GeV 



E > 25 GeV 




0.1 



0.2 



0.3 



Cluster Mass (GeV/c^ 



Figure 17: Comparison of the 7r° mass distributions of unresolved 7r°'s between data 
(points) and simulation (histogram) in different energy ranges. 
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• bremsstrahlung: bremsstrahlung photons radiated along the final charged particles 
in T decay, including the contribution from the detector material for electrons, 

• radiative: initial and final state radiation non-coUinear to charged particles in the 
final state, 

• genuine photons from 7r° decays where the partner photon is lost because of energy 
threshold, reconstruction inefficiency, cracks or overlap with another electromagnetic 
or hadronic shower, 

• genuine single photons, mostly from cu — > 7r°7 and 77 — > 77, and 

• fake photons. 

With respect to the previous analysis, the cut on Pj is tightened for residual 7r°'s since 
better discrimination between real and fake photons is now achieved. The cut used, 
InP^citt > —0.23 — 0.40 x Ey, has a rapidly increasing efficiency above 50% below 0.5 GeV 

and is about 90% at 3 GcV. 

The radiative and bremsstrahlung photons arc selected using the same method as 
described in the previous analysis [3]. Estimators Pbrem, Prad and P^o^^ are calculated to 
select photons from bremsstrahlung, radiative processes and from tt*^ decays, respectively. 
To compute these estimators, the angle between the photon and the most energetic 
charged track is used, in addition to the discriminating variables considered previously 
for the photon identification. Radiative and bremsstrahlung photons are not used in the 
r decay classification discussed below. 

The behaviour of these estimators was studied in Ref . [3] . The agreement between the 
number of bremsstrahlung photons in data and simulation is good, and affects mainly the 
electron r decay channel where this contribution is important (however it does not affect 
the rate). Bremsstrahlung photons in hadronic channels (i.e. radiation along the final 
state charged pions) are at a much lower level and are difficult to pick up unambiguously 
in the data. The estimate of the effect of this contribution largely relies on the description 
of radiation at the generator level in the Monte Carlo [22]. This point is addressed in 
Section 10.11. 

7 Decay classification 

Each r decay is classified topologically according to the number of charged hadrons, the 
charged particle identification and the number of 7r°'s reconstructed. While for one-prong 
and five-prong channels the exact multiplicity is required, the track number in three-prong 
channels is allowed to be 2, 3 or 4, in order to reduce systematic effects due to tracking 
and secondary interactions. Thus 13 classes are defined as given in Table 3. In this 
table, the right-most column shows how the different considered r decays contribute to 
the reconstructed channels as defined. Throughout this paper only t~ decays are cited 
and charge-conjugate modes are implied. 
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Table 3: Definition of the reconstructed r decay classes. All r decay modes implemented 
in the simulation are specified for each class. The notation r stands for r~ and the charge 
conjugate states are implied, while h stands for any charged hadron (tt or K). 



11 1 

Class label 


Reconstruction 
criteria 


Generated r decay 


£3 












h 


1 /i 


T ^ K Ur n 
T ^ K Ut- 




1 /i + 7r° 


T — > p T — > K TT K Uj- 

T ■K'n^'K t ^ K*~ v-j. 


h 27r° 


1 /i + 27r° 


- - i2\ 

T — > a^^ Vt t — > to Vt 

/ ' -£1- -^/l I^T 1 ' -li I^T 


h 37r° 


1 /i + 37r° 


T 7r-37r° T ^ K-n^K° z/^ 


h 47r° 


1 /i + > 47r° 


T ^ 7r-47r° i/^ _ 


3/i 


2 -Ah 


T ^ a{ Vt t — > K~K^'K~ Vt 

K*~ Ur T ^ n-K^K° Ur 
T K~7l^n~ Ur T — > K~K^ Ut 


3/i 7r° 


2 - 4/i + 7r° 


r — >• TT TT ii Vt- 


3/i 27r° 


3/i + 27r° 


T 27r-7r+27r" i/^ ^"^^ „ 

T — > TT A^A i/,- 


3/i 37r° 


3/i + > 37r° 


r ^ 27r'~7r"'"37r° i/^ 


5/i 


5/i 


T — > 37r~27r"'" i/^ r ^ 7r~i^°X° z/,- 


5/i 7r° 


5/i + 7r° 


T — > 37r^27r+7r° z/^- 



2 With a; ^ 77% 
^ With 77 — > 77 
^ With 77 ^ SttO 
^ This channel includes r 
^ This channel includes r 
With T] — > 7r~7r+7 
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TTW Vt with W — > TT TT+TT*^ 

TTTT^o; Vr with w — > 7r~7r+7r' 



The definition of tlie leptonic cliannels requires an identified electron or muon witli 
any number of photons. Some cuts on the total final state invariant mass are introduced 
to reduce feedthrough from hadronic modes [2]. Also some decays with at least two 
good electron tracks are now included, when one or more of such tracks are classified as 
converted photons. 

In the previous analysis [3], the 3/i27r°, ShSir^ and hAir^ channels (where h stands for 
any charged hadron: tt or K) suffered from large backgrounds and consequently had a 
low signal-to-noise ratio. Most of these backgrounds are due to secondary interactions 
of the hadronic track with material in the inner detector part. In order to improve the 
definition of these channels the following steps are taken: (1) require the exact charged 
multiplicity rich = 3 (instead of 2, 3, or 4) for 3/i27r° and ShSn^, (2) demand a maximum 
impact parameter of charged tracks less than 0.2 cm (instead of 2 cm) for ShSn^, and 
(3) require the number of resolved 7r°s in 3/i37r° to be 3 or 2, and 4 or 3 in Mtt^. With 
these tightened cuts the signal-to-noise ratio improves significantly with a small loss of 
efficiency. 

It should be emphasized that all hemispheres from the selected rr event sample are 
classified, except for single tracks going into an ECAL crack (but for identified muons) 
or with a momentum less than 2 GeV (except for identified electrons and hemispheres 
with at least one reconstructed tt"). These latter decays, in addition to the rejected 
ones in the 3/i27r° or 3/i37r° channels are put in a special class, labelled 14, which then 
collects all the non-selected hemispheres. By definition, the sample in class-14 does not 
correspond to one physical r decay mode. In fact, it follows from the simulation that this 
class comprises about 21% electron, 27% muon, 41% 1-prong hadronic and 11% 3-prong 
T decays. However, consideration of class-14 events can test if the rejected fraction is 
correctly understood, as discussed later. 

The numbers of r's classified in each of the considered decay channels are listed in 
Table 4. 

The KORALZ07 generator [17] in the Monte Carlo simulation incorporates all the 
decay modes considered in Table 3. Since the Mtt^ decay channel is not included in the 

standard Monte Carlo, a separate generation was done where one of the produced r is 
made to decay into that mode using a phase space model for the hadronic final state, while 
the other r is treated with the standard decay library. The complete behaviour between 
the generated decays and their reconstructed counterparts using the decay classification 
is embodied in the efficiency matrix. This matrix Sj^ gives the probability of a r decay 
generated in class j to be reconstructed in class i. Obtained initially using the simulated 
samples, it is corrected for effects where data and simulation can possibly differ, such as 
particle identification, as discussed previously, and photon identification as affected by 
the presence of fake photons. The latter correction is presented next, before exhibiting 
the final corrected matrix. 
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Table 4: Number of reconstructed events in 1991-1993 and 1994-1995 data sets in the 
different considered topologies. 



Reconstr. class 


nf (91-93) 


nf' (94-95) 


e 


22405 


33100 




22235 


32145 


h 


15126 


22429 




32282 


49008 




12907 


18317 


/i37rO 


2681 


3411 




458 


499 


3h 


11610 


17315 


3/i7r° 


6467 


9734 




1091 


1460 


3/i37r° 


124 


150 


5h 


60 


105 


5/i7r° 


36 


59 


Class 14 


4834 


7100 


sum 


132316 


194832 



8 Adjusting the number of fake photons in the 
simulation 

The number of fake photons in the simulation does not agree well with the rate in data. 
Such a discrepancy is seen in the comparison between the photon energy distributions in 
data and simulation given in Fig. 18: at low energy there is clear indication that the rate 
of simulated fake photons is insufficient. This effect will affect the classification of the 
reconstructed final states in the Monte Carlo and bias the efficiency matrix constructed 
from this sample. A procedure must be developed to correct for this effect. 

8.1 The method 

Taking explicitly into account the number k of fake photons in a r decay, the efficiency 
matrix can be rewritten as 

^ji — '^^jik • ^jk (6) 

k 

where ejik is the efficiency for a produced event in class j with k fake photons to be 
reconstructed in class i. It can be calculated from the simulation directly, given the 
number of fake photon in each event. Above ujjk is the fraction of produced class j events 
with k fake photons. The insufficient fake photons in simulation is due to the fact that 
ufl'^ is different from uff". 
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Figure 18: Energy spectra of observed photons in data (error bars) and Monte Carlo 
(histogram) samples. When normalized above 5 GeV photon energy the two distributions 
show a clear discrepancy in the low energy region. The shaded histogram shows the 
contribution from non-r backgrounds. 



Assuming each fake photon is produced randomly, one can get the fake photon 
multiplicity in Monte Carlo by in principle randomly removing a fraction of / fake photons 
in real data sample, i. e. 

uo^^ = Acu'^"*" (7) 

where a;''"*" = {ujo,uji,uj2, ■ ■ ■ ,ujn)data and u^"^ = {ujo,uji,uj2, ■ ■ ■ ,uJn)MC are the fake 
photon multiplicity distributions in data and simulation respectively. The matrix A is 
given by 

Am = Ct'f-'{1 - ff (8) 

if < /, otherwise it is 0. Here /c, I run from up to the maximum number of fake 
photons. The above equation holds for all the produced classes, with one parameter /, 
for each class. 

The number of fake photons in data can be expressed as 

■\TfAata T^Tproduced,data i V^/ /i— 1\ MC /r\\ 

K =2^NJ 2^£jik - k-}^{Aj )h-^jz^ (9) 

j k I 

i is the index for reconstructed class. In this set of equations, the numbers iv/''^"*'^ are 
obtained from photon probability fits in the reconstructed data samples, f^Y"'^^'^^'^''''"'^"' 
from the produced rr number and the measured branching ratios (implying an iterative 
procedure which in fact converges extremely fast), Sjik from the simulation with a proper 
procedure to determine k (called matching in the following), ujj{'" from the simulation 



35 



and matching, while {Aj depends on the fj unknowns to be determined by solving 
the equations. 

Thus to measure the fj factors, the number of fake photons in each simulated event 
and the number of fake photons in each reconstructed class in data must be determined. 



8.2 Counting fake photons in simulation: the photon matching 
procedure 

To count the number of fake photons in each simulated event, the philosophy is to remove 
all the clusters reconstructed in ECAL associated with produced photons at the generator 
level. The remaining clusters are then declared to be fake photons. 

To match a produced photon to an ECAL cluster, the cluster with minimum is 
selected, with 




where the quantities with S correspond to differences between reconstructed and generated 
levels, while the cr denote the corresponding expected resolutions. A matching flag M 
is assigned to each cluster according to the value. The photon sample with M = 1 
corresponds to small should be almost exclusively composed of genuine photons, 

with a small feedthrough from fake photons. The sample with intermediate is still 
dominated by genuine photons with a larger fake photon content. At the end of the loop 
over the produced (genuine) photons the unmatched clusters are given a flag M = — 1, 
corresponding to an almost pure sample of fake photons with a small amount of genuine 
photons. 

Many effects affect the matching procedure, leading to an incorrect estimate of the 
number of fake photons in the event: energy threshold, ECAL non-sensitive regions 
(cracks, acceptance), merging with other neutral clusters, merging with charged tracks, 
accidental matching. These effects have been studied very carefully and cuts are chosen 
so that the real photon contamination in the fake photon sample is the same as the fake 
photon contamination in the real photon sample. In this way the number of fake photons 
is given by the number of M = —1 photons, facilitating considerably further treatment 
of the fake photon contribution. 



8.3 Fits to data 

Using as reference the P-y distributions for M = — 1 and M 7^ —1 of a reconstructed class 
in the simulation, the distribution of all photons in a class i of data is fitted with 
these two components, and then corrected by the fraction of contaminated parts in the 
two distributions. Figure 19 shows a typical fit in this procedure. The fitted fake photon 
numbers in each observed class after correction are listed in Table 5, with the error from 
the fit and the correction. 
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Figure 19: Fit to the photon probabihty distribution in observed /ivr*^ events for the 1994- 
1995 data set. From top to bottom, are the probability distributions of fake and real 
photon samples in the simulation from the matching procedure, and the fit to the data 
distribution using the two Monte Carlo components. 
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Table 5: Number of fake photons in 1991-93 and 1994-1995 data sets. 



Reconstr. 


class 91-93 


94-95 


e 


2448 ± 75 


3552 ± 110 


1^ 


183 ± 24 


170 ± 26 


h 


3718 ± 58 


5766 ± 98 




9708 ± 196 


14844 ± 350 




5941 ± 172 


8837 ± 237 


/i37r° 


2311 ± 83 


3085 ± 90 




574 ± 35 


677 ± 35 


3h 


5769 ± 88 


8639 ± 162 


3/i7rO 


4130 ± 98 


6444 ± 112 


3/i27r° 


1129 ± 48 


1432 ± 52 


3/i37r° 


183 ± 18 


225 ± 24 


5h 


45 ± 6 


73 ± 13 




30 ± 6 


52 ± 8 


Class 14 


1369 ± 43 


2014 ± 62 



8.4 An independent check 

The total excess of fake photons in data with respect to Monte Carlo from the fits to 
the distributions in each observed class can be checked comparing the photon energy 
spectra of data and Monte Carlo, after normalization to the number of r events and 
correcting the simulation for measured branching ratios. The latter method is justified 
by the fact that fake photons show up as an excess at energies typically less than 1-2 
GeV. The former estimate yields 7279±334 and 11419±511 for 1991-1993 and 1994-1995 
data samples respectively summing up all channels, while the latter yields 7007 ±307 and 
10820 ± 376 for the two data sets. The differences 

A91-93 = 272 ±454 (11) 
A94_95 = 599 ±634 (12) 

indicate that the number of fake photons estimated from the fit procedure is reliable. 

8.5 Solving for the fake photon correction factors 

In order to reduce the number of parameters to be determined, it is assumed that the 
correction factors of some channels with small branching ratios are the same as those of 
the nearby channels with comparable multiplicity. Also the /i channel is not considered 
since it is hardly influenced by fake photons and the feedthrough to other channels is very 
small. Therefore the correction factors are only fitted for e, h, hir^, h2TT'^, 3h and 3/t7r° 
channels. The results are given in Table 6. They are consistent for the two data/Monte 
Carlo sets. The resulting xV degrees of freedom (DF) is 12.4/5 and 15.2/5 for 1991-1993 
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Table 6: Fake photon correction factors for 1991-1993 and 1994-1995 data sets. 



produced class 


91-93 (%) 


94-95 (%) 


e 


41.0 ± 1.9 


36.1 ± 2.2 


h 


8.3 ± 2.2 


12.1 ± 2.3 




21.7 ± 2.2 


20.5 ± 2.6 




13.0 ± 2.7 


18.7 ± 2.0 


3h 


12.9 ± 1.6 


14.0 ± 1.9 


3/i7r° 


25.1 ± 2.5 


27.4 ± 1.6 



and 1994-1995 data respectively. The final error on the fitted correction factors has been 
properly enlarged to take into account these somewhat large values. 



8.6 Comparison to the previously used method 

In the published analysis using 1991-1993 data [3] a much less sophisticated approach was 
taken. The deficit of fake photons in the simulation was determined in a global way to be 
(16 ± 8)%, common to all channels. Since it would have been dehcate to generate extra 
fake photons, the procedure chosen was to actually do the opposite, i.e. randomly kill 
identified (in the sense of the matching procedure discussed above) fake photons in the 
simulation, determine the new efficiency matrix and use the deviations to correct in the 
opposite direction. 

The current way of dealing with the fake photon problem is both more precise and 
more reliable. The fact that different channels are treated separately provides a handle 
on the different origins of the fake photons, since, for example, fake photons in the 
h class only originate from hadronic interactions, whereas they come from both hadronic 
interactions and photon shower fluctuations in the hn^ class. Also the previous procedure 
was dependent on the quality of the photon matching in simulation and the killing of fake 
photons entailed the possibility to mistakenly removing some genuine photons as well. 



9 Determination of the branching ratios 

The branching ratios are determined using 

nf-nf^ = Y.SjiNp" (13) 
j 

njprod 

Bi = 7 (14) 

where nf' is the number of observed r candidates in reconstructed class i, n^^ the 
non-T background in reconstructed class i, Sji the efficiency of a produced class j event 
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Table 7: Efficiency matrix for 1994-95 data (in %). Generated classes are given in the 
first row, and reconstructed classes in the first column. All the corrections are applied 
except that for nuclear interactions and hadron misidentification in multiprong channels. 
The last row gives the selection efficiency for each produced class. 





e 


M 


h 




h2^'> 




h47r" 


3h 


3hn'> 


3h27r" 


3;i37r" 


5fe 


5)1,^" 


e 


73.26 


0.01 


0.41 


0.45 


0.34 


0.25 


0.74 


0.02 


0.02 


0.05 


0.00 


0.00 


0.00 




0.01 


74.49 


0.63 


0.22 


0.07 


0.21 


0.33 


0.01 


0.01 


0.00 


0.00 


0.00 


0.00 


h 


0.26 


0.76 


65.03 


3.66 


0.34 


0.06 


0.00 


1.44 


0.10 


0.08 


0.00 


0.80 


0.00 




1.02 


0.26 


4.70 


68.19 


11.31 


2.16 


0.49 


0.48 


1.28 


0.62 


0.05 


0.24 


0.00 


h2-n-° 


0.12 


0.01 


0.33 


5.67 


57.68 


23.13 


7.57 


0.08 


0.39 


1.48 


0.24 


0.04 


0.00 


Mir" 


0.01 


0.00 


0.07 


0.41 


6.92 


43.06 


38.15 


0.01 


0.10 


0.37 


0.71 


0.04 


0.00 


h4TT° 


0.00 


0.00 


0.02 


0.05 


0.67 


6.25 


25.26 


0.00 


0.02 


0.11 


0.19 


0.00 


0.00 


Sh 


0.01 


0.02 


0.2,-) 


0.07 


0.03 


0.00 


0.00 


67.98 


6.77 


0.80 


0.03 


22.11 


2.52 


Shir" 


0.01 


0.01 


0.22 


0.56 


0.27 


0.06 


0.06 


7.29 


58.90 


16.53 


4.46 


7.07 


16.04 


3h27i-° 


0.00 


0.00 


0.04 


0.06 


0.10 


0.08 


0.02 


0.41 


6.02 


40.42 


25.02 


0.28 


0.65 


3h3n° 


0.00 


0.00 


0.00 


0.00 


0.00 


0.01 


0.05 


0.02 


0.41 


6.19 


28.98 


0.00 


0.00 


5h 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.01 


0.01 


0.00 


0.00 


38.70 


4.58 




0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.01 


0.02 


0.03 


0.08 


2.99 


38.72 


Class 14 


3.27 


4.17 


6.38 


0.73 


1.08 


1.71 


1.75 


0.80 


3.66 


9.96 


13.87 


5.03 


9.75 


sum 


77.06 


79.72 


78.08 


79.97 


78.81 


76.97 


74.42 


78.66 


77.71 


76.64 


73.64 


77.30 


72.26 



reconstructed as class i, and Nj^" the produced events number of class j. The class 
numbers i, j run from 1 to 14, the last one corresponding to the rejected r candidates. 

The efficiency matrix Sji is determined from the Monte Carlo, but corrected using 
data for many effects such as particle identification efficiency, rr selection efficiency and 
fake photon simulation. It is given in Table 7 where the last line shows that the selection 
efficiency remains constant within ± 5% for all 13 considered topologies. 

It should be noted that the efficiency matrix eji is independent of the r branching 

ratios used in the simulation, except for the subclasses contributing to each defined class 
as shown in Table 3. The effect depends however on small branching ratios for final states 
including kaons and the procedure used for this correction relies on the ALEPH measured 
values [9]. 

The analysis assumes a standard r decay description. One could imagine unknown 
decay modes not included in the simulation, but since large detection efficiencies are 
achieved in the rr selection which is therefore robust, so that these decays would be 
difficult to pass unnoticed. An independent measurement of the branching ratio for 
undetected decay modes, using a direct search with a one-sided r tag, was done in 
ALEPH [23], limiting this branching ratio to less than 0.11% at 95% CL. This result 
justifies the assumption that the sum of the branching ratios for visible r decays is equal 
to unity. 

The equations are conveniently solved by a minimization technique and all corrections 
not yet included in the efficiency matrix are applied. These small final corrections are 
discussed in the following section on systematic studies. The branching ratios are obtained 
and listed in Table 8. The results on the class- 14 "branching ratio" are discussed in 
Section 13.1. 
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Tabic 8: Branching ratios (%) from 1991-1993 and 1994-1995 data sets; the first error is 
statistical and the second is systematic. 



Topology 


91-93 


94-95 


e 


17.859 ± 0.112 ± 0.058 


17.799 ± 0.093 ± 0.045 




17.356 ± 0.107 ± 0.055 


17.273 ± 0.087 ± 0.039 


h 


12.238 ± 0.105 ± 0.104 


12.058 ± 0.088 ± 0.083 




26.132 ± 0.150 ± 0.104 


26.325 ± 0.123 ± 0.090 


/i27rO 


9.680 ± 0.139 ± 0.124 


9.663 ± 0.107 ± 0.105 




1.128 ± 0.110 ± 0.086 


1.229 ± 0.089 ± 0.068 


MvrO 


0.227 ± 0.056 ± 0.047 


0.163 ± 0.050 ± 0.040 


?>h 


9.931 ± 0.097 ± 0.072 


9.769 ± 0.080 ± 0.059 




4.777 ± 0.093 ± 0.074 


4.965 ± 0.077 ± 0.066 


3/i27r° 


0.517 ± 0.063 ± 0.050 


0.551 ± 0.050 ± 0.038 


3/i37r° 


0.016 ± 0.029 ± 0.020 


-0.021 ± 0.023 ± 0.019 


5/i 


0.098 ± 0.014 ± 0.006 


0.098 ± 0.011 ± 0.004 




0.022 ± 0.010 ± 0.009 


0.028 ± 0.008 ± 0.007 


Class 14 


0.017 ± 0.043 ± 0.042 


0.099 ± 0.035 ± 0.037 



10 Determination of systematic uncertainties 
10.1 Methodology 

Wherever possible the efficiencies relevant to the analysis have been determined using 
ALEPH data, either directly on the rr sample itself or on specifically selected control 
samples, as for example in the case of particle identification. The resulting efficiencies are 
thus measured with known statistical errors. 

In some cases the procedure is less straightforward and involves a model for the 
systematic effect to be evaluated. An important example is given by the systematics 
in the simulation of fake photons in ECAL. In such cases the evaluation of the systematic 
error not only takes into account the statistical aspect, but also some estimate of the 
systematics involved in the assumed model. The latter is obtained from studies where the 
relevant parameters are varied in a range consistent with the comparison between data 
and Monte Carlo distributions. 

Quite often a specific cut on a given variable can be directly studied. The comparison 
between the corresponding distributions, respectively in data and Monte Carlo, provides 
an estimate of a possible discrepancy whose effect would be to change the assumed 
efficiency of the cut. If a significant deviation is observed, a correction is applied to 
the simulation to obtain the nominal branching ratio results, while the error on the 
deviation provides the input to the evaluation of the systematic uncertainty. The analysis 
is therefore repeated with a full re-classification of all the measured r decay candidates, 
changing the incriminated cut by one standard deviation. Since the new samples slightly 
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Tabic 9: Selection efficiencies (%) for rr events, corresponding to cuts 5-9 described in 
Section 3.1 , measured by the break-mix method for the 1994-1995 data. 



class 


data (BM) 


Monte Carlo (BM) 


R=data/MC 


e 


95.70 ± 0.10 


95.54 ± 0.07 


1.0017 ± 0.0012 




96.10 ± 0.10 


96.20 ± 0.06 


0.9990 ± 0.0012 


h 


96.03 ± 0.10 


95.95 ± 0.07 


1.0009 ± 0.0013 




97.07 ± 0.10 


97.26 ± 0.04 


0.9980 ± 0.0011 


h2n^ 


95.36 ± 0.16 


95.68 ± 0.10 


0.9967 ± 0.0020 


3h 


94.77 ± 0.18 


95.10 ± 0.10 


0.9965 ± 0.0022 


3hn° 


92.53 ± 0.22 


92.93 ± 0.19 


0.9957 ± 0.0031 



differ from the nominal ones because of feedthrough between the different channels, the 
modified results are affected both by the systematic change in the variable value and the 
statistical fluctuation from the event sample which is uncommon to both samples. In this 
case the final systematic uncertainty is obtained by adding linearly the modulus of the 
systematic deviation observed and the statistical error from the monitored independent 
sample. This procedure is followed in a systematic way for the studies listed below and 
its description will of course not be repeated each time. 

Finally, the systematic deviations for each study are obtained with their sign in each 
measured decay channel, thus providing the full information on the correlations between 
the results and allowing the corresponding covariance matrix to be constructed. 

10.2 Selection of rr events 

Selection efficiencies of all the dominant channels have been determined from data using 
the break-mix method as discussed in Section 3.2. The efficiencies are very similar for the 
different data sets. They are listed in Table 9 for the 1994-1995 sample. 

The obtained values are then allowed to float around the nominal values with their 
corresponding errors, and 1000 toy experiments are generated to estimate the effect on 
the branching ratios. The errors are listed in Table 11 in the column 'sel'. 

10.3 Non-T background 

Most non-T background contributions in each r decay channel have been measured directly 
in data following the procedure described in Section 3.3. The contamination from hadronic 

Z decays has been determined by Monte Carlo simulation. In that case the relative error 
of non-r background is taken as 30%. Toy experiments are used where new values for the 
non-r background are generated according to Gaussian distributions with the measured 
error. Branching ratios are reevaluated in each case and the resulting standard deviations 
are quoted as errors. They are given in Table 11 in the column 'bkg'. 
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10.4 Particle identification 



The same technique is used to estimate the systematic error on particle identification and 

misidentification efficiencies which are obtained from measurements on control samples 
(sec Section 4.2) with their own statistical and systematic uncertainties. The errors on 
the branching ratios are listed in Table 11 in the column 'pid'. 

10.5 Photon detection efficiency 

Photon detection efficiency at low energy and high energy are studied by different means. 

10.5.1 Photon efficiency at low energy 

The number of real photons in the data sample is obtained from the fit of the 
distribution in each energy bin with the distributions of genuine and fake photons from 
Monte Carlo sample. The energy spectra of genuine photons for both data and Monte 
Carlo samples are then compared in order to obtain the relative photon efficiency between 
data and Monte Carlo (Fig. 20). Some discrepancy is observed at low energy, indicating 
that the efficiency in the data sample is different from the simulation. It should be 
emphasized that this effect originates almost exclusively from the detector simulation 
which was improved for the 1994-1995 sample. It is taken into account in the analysis 
by randomly killing Monte Carlo photons in the corresponding energy bins. The error on 
the efficiency is propagated to the final branching ratios and are given in Table 10 in the 
column 'eff'. 

One could question whether the above effect is indeed related to the simulation of 
photon efficiency or rather induced by different physics in data and in simulation. The fact 
that the first explanation is the correct one is made clear by the following observations: (1) 
the discrepancy occurs only near photon-detection threshold with a characteristic shape, 
(2) the branching ratio values are adjusted in the Monte Carlo to match those measured 
in data, (3) the effect is observed in all photonic classes, and (4) the effect is quantitatively 
different in the barrel and endcap parts of ECAL. 

The efficiency is measured for 1991-1993 and 1994-1995 samples separately, and the 
values are found to be different for 1-prong and 3-prong r decays, which is reasonable 
considering the different environment of hadronic interactions in ECAL. The efficiencies 
in the ECAL barrel and endcaps are also measured separately, because the photon 
performance is quite different in the endcaps due to more material present in front of 
the calorimeter. 

10.5.2 Photon efficiency at high energy 

A loss of high energy photons occurs because of the merging of their showers into nearby 
charged track clusters. Photon reconstruction requires a cut on dc at 2 cm in order to 
separate the photon candidate from the charged track cluster. The corresponding loss 
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Figure 20: Relative photon efficiency in data and Monte Carlo for 1991-1993 (open circles) 
and 1994-1995 (black dots) data sets. The plots are separated for the barrel and endcap 
sections of the calorimeter, and for one- and three-prong r decays. 



Table 10: Systematic error for photon and tt" reconstruction in 1994-1995 data sample. 
All numbers are absolute branching ratios in %. The labels are defined as follows: photon 
efficiency at low energy (eff), photon efficiency at high energy (dgt), converted photons 
(cnv), photon identification efficiency (prb), photon energy calibration (cal), fake photon 
correction (fak), 7r° reconstruction efficiency (dij), bremsstrahlung and radiative photons 
(bms). The total 7r° systematic uncertainty is given in the rightmost column. 



Topology 


eff 


dgt 


cnv 


prb 


cal 


fak 


dij 


bms 


total TT 


e 


0.006 


0.002 


0.005 


0.000 


0.004 


0.003 


0.005 


0.003 


0.011 


M 


0.002 


0.001 


0.001 


0.000 


0.002 


0.000 


0.001 


0.002 


0.004 


h 


0.022 


0.029 


0.006 


0.056 


0.009 


0.011 


0.002 


0.019 


0.071 




0.015 


0.011 


0.009 


0.024 


0.011 


0.048 


0.006 


0.023 


0.063 


h2iT° 


0.033 


0.013 


0.013 


0.041 


0.018 


0.048 


0.025 


0.038 


0.089 


/i37r° 


0.012 


0.019 


0.010 


0.035 


0.014 


0.008 


0.030 


0.013 


0.056 


/i47r° 


0.017 


0.012 


0.005 


0.012 


0.006 


0.010 


0.007 


0.007 


0.029 


3h 


0.028 


0.012 


0.007 


0.029 


0.010 


0.012 


0.005 


0.011 


0.046 


3/177° 


0.004 


0.003 


0.010 


0.014 


0.006 


0.020 


0.010 


0.017 


0.033 


3/i27r0 


0.016 


0.010 


0.008 


0.012 


0.006 


0.008 


0.003 


0.007 


0.027 


3/i37r° 


0.002 


0.004 


0.003 


0.004 


0.003 


0.002 


0.005 


0.004 


0.010 


5h 


0.001 


0.001 


0.000 


0.001 


0.000 


0.000 


0.000 


0.001 


0.002 




0.001 


0.000 


0.001 


0.001 


0.001 


0.000 


0.001 


0.001 


0.002 


Class 14 


0.008 


0.004 


0.004 


0.002 


0.003 


0.002 


0.004 


0.006 


0.013 
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Figure 21: Comparison of d^. distributions for high energy photons [E^ > 3 GeV) in data 
(points) and in Monte Carlo (histogram) for 1991-1993 and 1994-1995 data sets. The 
distributions are normahzed to each other for > 8 cm. 

of efficiency depends on the simulation of hadronic interactions in ECAL, which can 
be checked against data by comparing the actual d^ distributions. Figure 21 shows the 
comparison for photons with energy greater than 3 GeV. After normalization to dc> 8 cm, 
the fraction of photons below 8 cm is (0.30 ± 0.15)% more in data for 1991-1993 and is 
(0.37 ± 0.13)% more in data for 1994-1995. 

The observed difference is applied in the analysis and varying it by one standard 
deviation provides an estimate of the systematic error induced by the dc cut. The errors 
are listed in Table 10 in the column 'dgt'. 

10.6 Converted photons 

The rate of reconstructed converted photons depends on the conversion probability and the 
misidentification of hadronic tracks in a multihadron environment. The misidentification 
of protons originating from nuclear interactions should also be considered in this context 
since their energy loss by ionization can be close to that expected for electrons in the low 
momentum range. 
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Figure 22: Ratio of converted photon energy spectra in data and Monte Carlo. 



Figure 22 shows the ratio of converted photon spectra in data and Monte Carlo, 
indicating an inadequate simulation of the rate of converted photons. Since in the analysis 
converted photons are given an identification probability of 1, the net effect is an efficiency 
increase for photons in data. To correct for this, some converted photons in the data 
sample are randomly declared to be calorimeter photons according to the observed excess 
and given an identification probability value according to the corresponding ECAL photon 
distribution. Their fate in the analysis thus follows that of the ordinary calorimetric 
photons which can be removed when probability cuts are applied. 

The error on the conversion rate ratio is propagated to branching fractions through a 
full reanalysis. The corresponding systematic errors are listed in Table 10 in the column 
'cnv'. 

The misidentification of hadrons to electrons in a multihadron environment is 
measured by comparing events reconstructed as 3/i7r° with two tracks identified as hadrons 
and one reconstructed converted photon: this sample is dominated by produced 3h events 
with one hadronic track misidentified as electron. Comparing the event numbers in data 
and Monte Carlo after normalization to the total r event number, the rate of hadron 
misidentification is found to be (0.06 ±0.03)% lower in simulation for the 1994-1995 data 
sample, corresponding to an absolute change in the 3h fraction of (+0.010 ± 0.005)%. A 
bigger difference is observed in the 1991-1993 data sample where the misidentification 
probability is found to be (0.25 ± 0.04)% lower in the Monte Carlo. This problem 
only affects feedthrough between nearby 3-prong channels, the largest correction being 
to the 3h channel, (+0.034 ± 0.007)%, while the effect on 3/i7r° and 3/i27r° modes is 
(-0.019 + 0.004)% and (-0.015 + 0.003)%, respectively The found difference is illustrated 
in Fig. 23 where a comparison of the data and Monte Carlo distributions of the momenta 
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Figure 23: Comparison of hadron misidentification in multihadronic final states (see 
text for details). Data are given by points and Monte Carlo by histogram, while the 
shaded histogram shows the contribution from non-r background. A clear improvement 
is observed in the 1994-1995 simulation. 

of the converted photons is shown. The error of this correction is included under 'cnv' in 
Table 10. 

Proton tracks are observed when hadrons interact with the material of the inner 
detectors. Since at low energy, only dE/dx information is used for particle identification, 
protons can be misidentified as electrons for momenta between 0.9 and 1.1 GeV, as seen 
in Fig. 24. By selecting tracks with large dE/dx at even lower momenta so that protons 
dominate the sample, the proton rates can be compared in data and in simulation directly. 
Good agreement is observed and the systematic error on r branching ratios from proton 
tracks in the photon conversions is negligible. 

10.7 Photon identification efficiency 

Since the reference distributions for photon identification are not exactly the same for 
data and Monte Carlo (even after corrections), the same probability cut will have different 
efficiencies for the two cases. Figure 25 shows the comparison of distributions between 
data and simulation for residual photons in the 0.3-3.0 GeV range: good agreement is 
observed overall, however some small shape differences can be seen. A comparison of 
the same distributions is performed for photons in small energy slices, after correcting 
the branching ratios in the Monte Carlo and for the non-r background. The numbers of 
photons surviving the cut are compared, and the efficiency is found to be a little larger in 
the simulation. This effect is then corrected in the analysis by using a slightly tighter cut 
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Figure 24: Measured dE/dx as a function of momentum for all (good and bad) tracks in 
the selected r sample: the proton band from nuclear interactions is clearly observed, and 
predicted by the simulation. 
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Figure 25: Comparison of photon identification probability distributions in data (points) 
and Monte Carlo (solid histogram, after fake photon correction) for residual photons in 
the 0.3-3.0 GeV energy range. The dashed histogram shows the simulated distribution 
before fake photon correction, and the shaded histogram indicates the contribution of 
non-r background. 



for the simulation, corresponding to the observed difference in the numbers of accepted 
photons. The systematic errors on branching ratios appear in Table 10 in the column 
'prb'. 



10.8 Photon energy calibration 

The photon energy scale relative to the simulation is calibrated using high energy electrons 
from r decays as discussed in Section 5.2.5. The parameterization of the E/p ratio is 
extrapolated to low energy for both data and Monte Carlo. The error on the calibration 
factors is propagated to branching ratios and are listed in Table lO('car). 
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10.9 Fake photon correction 

The error on the fake photon correction comes from the uncertainty in measuring the 
total number of fake photons in each observed class for data, and in determining the fake 
photon multiplicity in each generated Monte Carlo event, using the matching procedure 
described in Section 8.2. The latter has a negligible effect on branching ratios since 
the dependence of the efficiency matrix on the fake photon multiplicity is weak. The 
fake photon numbers are determined from fits to the distributions, and the errors are 
obtained in the global fit of the fake photon correction factors. These errors are then 
propagated to branching ratios. To reduce the total number of parameters the correction 
factors on low branching ratio channels are assumed to be equal to those of the nearby 
high branching ratio channels, thus somewhat increasing the of the fits. The fitted 
errors are then enlarged by the ^J'x^ / DF values in each channel. The errors are reported 
in Table 10 in the column 'fak'. 

10.10 7r° reconstruction efficiency 

To estimate the uncertainty on the reconstruction of the resolved vr^'s, a comparison of the 
Dij (Section 6) distributions is shown in Fig. 26 after branching ratio, fake photon, and 
non-r background corrections are applied in the simulation. Good agreement is observed. 
The agreement on the reconstruction efficiencies of resolved 7r°'s in data and Monte Carlo 
is tested with a precision of 0.17% and 0.14% for the 1991-1993 and 1994-1995 samples 
respectively. They can be reproduced by changing the value of the Dij cut by 0.14 x 10~^ 
and 0.12 x 10""^ for 1991-1993 and 1994-1995 respectively, the nominal cut b cing set at 
1 X 10~^. The effect of changing the cut by this amount in the simulation was taken as 
the systematic error, listed in Table 10 in the column 'dij'. 

10.11 Simulation of bremsstrahlung and radiative photons 

Figure 27 shows the comparison of the energy spectra of selected radiative and 

bremsstrahlung photons in data and simulation after normalization to the total number 
of r events. Reasonable agreement is found except for the radiative photon spectrum in 
the 1991-1993 sample, where more photons are found in the simulation. A comparison 
of 1991-1993 and 1994-1995 simulated samples clearly shows that too many radiative 
photons were produced in 1991-1993 at the generator level (Fig. 28). To correct for this, 
photons matched as radiative in the 1991-1993 Monte Carlo sample are deleted randomly 
according to the 1994-1995 distribution, which agrees with data. The angle between the 
radiative photon and the jet axis is also taken into account in this procedure. 

The statistical error in the measurement of the numbers of radiative and 
bremsstrahlung photons is propagated to the r branching ratios and listed in Table 10 in 
the column 'bms'. 

It turns out that most of the bremsstrahlung photons originate from radiation in the 
detector material by electrons in the decay t" — > i/^e~l7e. The comparison in Fig. 27 
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Figure 26: Comparison of Dij distributions in data (points) and simulation (dark 
histogram) for 1991-1993 and 1994-1995 data sets. The hght histograms show the 
distribution before fake photon correction. 
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Figure 27: Comparison of measured energy spectra of bremsstrahlung and radiative 
photons in data (points) and Monte Carlo (histogram) for 1991-1993 and 1994-1995 data 
sets. 
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Figure 28: Comparison of the angle between produced radiative photons and jet axis 
in 1991-1993 (points) and 1994-1995 (histogram) simulated samples for different energy 
slices. 
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shows that the description of radiation is adequate in the simulation. However, any 
difference at this level would only affect the electron momentum spectrum, but not the 
decay classification. 

The situation is more delicate for radiative photons in hadronic decay channels. From 
simulation, it is found that only 23% of these photons are recognized as such, while the 
purity in the identified radiative photon sample is about 55% (others photons are either 
fake or originating from 7r° decays). These values are weakly dependent on the decay 
channel. Since the simulation is relied upon for the unobserved radiative photons, it is 
important to validate the radiation procedure in the simulation. 

In KORALZ07 [17] a model for radiation based on PHOTOS [22] is implemented: 
charged particles independently emit photons according to a bremsstrahlung probability. 
A test of this model can be performed in the tttt^ channel where a complete calculation 
of radiative corrections exists [24]. If a photon is radiated by either the p or the charged 
TT, then the final state can be identified as tttt^ if the radiative photon is classified as a 
bremsstrahlung or radiative photon, or undetected. Otherwise such a decay would be 
(wrongly) classified in the 7r27r'^ class and the corresponding bias should be corrected for 
through the simulation. It is thus important to verify the radiation model against the full 
calculation. A good agreement is found for the radiative rate for photon energies E* above 
12 MeV in the r centre-of-mass frame, comparable to the detection threshold of 350 MeV 
in the e+e" frame: the simulation using PHOTOS yields a radiative branching ratio of 
(2.91 ± 0.04) X 10^^ to be compared to a value of 2.9 x 10~^ from the exact calculation. 
In both cases, the radiative yield includes photons emitted by the r lepton and by the 
charged p and tt in the decay. The simplified bremsstrahlung probability is however no 
longer adequate for hard radiation: the corresponding values for E* > 300 MeV are 
(0.90 ± 0.09) X 10~^ for the simulation and 1.8 x 10"^ for the calculation. Although 
significant, such a difference has no practical influence in this analysis, since it only 
affects a negligible part of the radiative photon spectrum. The comparison of the total 
radiative yield above detection threshold leads to an absolute systematic error of 0.005% 
on the branching ratio in the tttt'^ channel, a factor of 3 to 4 smaller than the estimated 
uncertainties quoted in Table 10 from the direct data/simulation comparison of the 
radiative photon spectra. The latter values are therefore retained as final uncertainties 
for this source. 



10.12 Simulation of tt^ Dalitz decays 



The simulation of the tt Dalitz decays is checked at the generator level. Good agreement 
is found with the world average branching ratio. The quality of the simulation of 7r° Dalitz 
decays can also be checked in the distribution of the reconstructed photon conversion point 
near interaction point which agree well with data (Fig. 16). The possible systematic effect 
in the simulation of tt*^ Dalitz decays on the r branching ratios is found to be negligible. 
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Tabic 11: Total systematic errors for branching ratios measured from the 1994-1995 
data sample. All numbers are absolute in per cent. The labels are defined as follows: 
photon and reconstruction (tt^), event selection efficiency (sel), non-r background 
(bkg), charged particle identification (pid), secondary interactions (int), tracking (trk), 
Monte Carlo dynamics (dyn), Monte Carlo statistics (mcs), total systematic uncertainty 
(total). 



Topology 




sel 


bkg 


pid 


int 


trk 


dyn 


mcs 


total 


e 


0.011 


0.021 


0.029 


0.019 


0.009 


0.000 


0.000 


0.015 


0.045 




0.004 


0.020 


0.020 


0.021 


0.008 


0.000 


0.000 


0.015 


0.039 


h 


0.071 


0.016 


0.010 


0.022 


0.022 


0.014 


0.000 


0.019 


0.083 


h7T° 


0.063 


0.027 


0.019 


0.011 


0.045 


0.009 


0.000 


0.027 


0.090 


h2TT° 


0.089 


0.021 


0.014 


0.004 


0.007 


0.003 


0.040 


0.028 


0.105 




0.056 


0.012 


0.015 


0.000 


0.008 


0.001 


0.008 


0.030 


0.068 


Mtt^ 


0.029 


0.005 


0.011 


0.000 


0.015 


0.000 


0.000 


0.019 


0.040 


5h 


0.047 


0.021 


0.018 


0.004 


0.012 


0.014 


0.006 


0.015 


0.059 


3hTr° 


0.033 


0.017 


0.029 


0.002 


0.041 


0.009 


0.007 


0.018 


0.066 


3/i27r° 


0.027 


0.008 


0.015 


0.000 


0.009 


0.003 


0.012 


0.014 


0.038 




0.010 


0.012 


0.002 


0.000 


0.002 


0.001 


0.010 


0.006 


0.019 


5h 


0.002 


0.000 


0.002 


0.000 


0.000 


0.001 


0.000 


0.003 


0.004 




0.002 


0.000 


0.006 


0.000 


0.000 


0.000 


0.000 


0.002 


0.007 


Class 14 


0.013 


0.003 


0.022 


0.002 


0.024 


0.000 


0.000 


0.011 


0.037 



10.13 Nuclear interactions 

Secondary interactions of charged hadrons produced from r decays with material in the 
inner part of the detector are studied with the help of an enriched sample. Interaction 
events are characterized by the presence of many bad tracks and proton tracks. By 
requiring more than three bad tracks in the hemisphere or at least one proton track 
identified with energy loss information, interaction events are selected with an efficiency 
of about 50% and a purity of 70%, as estimated from the simulation. 

Figure 29 shows the comparison of data to Monte Carlo of the observed r decays 
in each reconstructed class for this interaction-rich sample, after normalizing to the total 
number of r events. Globally it is found that the MC overestimates the rate for secondary 
interactions by about 10%, with some channel-to-channel variation. 

By comparing the total number of interaction events in data and simulation, the 
difference of interaction rates is measured and corrected, modifying the efficiency matrix 
slightly. The branching ratios are then reevaluated with this new efficiency matrix. The 
differences between the new branching ratios and those without correction are taken as the 
conservative estimate of the systematic error from the simulation of nuclear interactions. 
This procedure is justified by the fact that, although the rate differences can be well 
measured, it does not ensure that the interaction dynamics is properly simulated. As 
shown in Fig. 29, even after the apphed global rate correction the agreement between 
data and Monte Carlo, although much improved, remains not very satisfactory. 

Another approach to study the influence of secondary interactions relies on the 
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Figure 29: Comparison of the numbers of selected r decays with secondary interactions 
for each channel (see Table 3 in data (points) and Monte Carlo (histogram)). The shaded 
part gives the contribution from non-r background, mainly from Z ^ qq events. The 
secondary tracks found in leptonic channels (classes 1 and 2) originate from the opposite 
(hadronic) decay. 
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Figure 30: Comparison of the distributions of the maximum value for selected r decays 
with secondary interactions in data (points) and Monte Carlo (histogram), the shaded 
part shows the contribution from non-r background, mainly from Z —>■ qq events. The 
global interaction rate difference has been corrected for. 



distribution of the distance in the transverse plane between the extrapolated track and the 
interaction point. Figure 30 shows the corresponding comparison of data to simulation 
(in three-prong channels, the largest of the three distances is plotted). After the global 
interaction rate correction, good agreement is observed for both one-prong and three- 
prong channels for 1991-1993 and 1994-1995 data sets, thus corroborating the findings of 
the other study. 

The errors are listed in Table 11 in the column 'int'. 



10.14 Tracking 

The only situation where tracking efficiency is relevant occurs in a three-prong decay 
when two of the tracks are lost. This effect can be studied by looking at events with same 
sign hemispheres. In order to reduce the effect of secondary interactions and asymmetric 
photon conversions, the requirements of no proton track in both hemispheres, less than 
four bad tracks and no track classified as converted photon, are made. In this way 121 ±9 
events are obtained in the Monte Carlo sample after normalizing to the total r events 
in the 1994-1995 data sample. The efficiency of measuring same sign events resulting 
from tracking is found to be 55% in simulation with a purity of 60%, the remainder being 
dominated by one-prong r decays with nuclear interactions. The corresponding number 
from the data sample is 180 ± 13, significantly larger than the rate seen in the simulation. 
Taking into account the corrected simulation rate of nuclear interactions as determined 
in the preceding section, the data sample produces an estimated excess of (81 ± 24)% 
same sign hemispheres compared to the simulation. This value corresponds to a tracking 
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efficiency difference of (0.3 ± 0.1)% between data and simulation. Tliis correction is 
applied to the efficiency matrix element corresponding to produced three-prong r decays 
reconstructed as 1-prong, which increases by about 40%, only about half of the three- 
prong r decays reconstructed as one-prong are due to tracking efficiency, the another half 
originating from — * tt+tt" decays with their secondary vertex far from the interaction 
point. The efficiency of three-prong to three-prong decays is decreased accordingly. The 
error on the tracking efficiency is propagated to the r branching ratio measurement. 
The largest correction in the branching ratios due to tracking is in the h channel, which 
decreases by (0.047 ± 0.014)% absolutely, while the 3h channel correspondingly increases 
by the same amount. 

The same analysis gives (104 ± 7) same sign hemisphere events in 1991-1993 Monte 
Carlo sample, while (109 ± 10) events are observed in data, thus these are in good 
agreement. Following the same procedure, the excess of misreconstructed events in the 
data sample is (9±21)%. No correction is taken into account this time and the systematic 
uncertainty is evaluated in the same way as for 1994-1995 sample. 

The errors are shown in Table 11 in the column 'trk'. 



10.15 Dynamics 

Uncertainties in the dynamics of the hadronic r decays can lead to systematic biases when 
computing the efficiency matrix from the simulation. The fact that selection efficiencies 
are large and weakly dependent on the hadronic invariant mass is an important factor 
hmiting those biases. Nevertheless systematic checks have been performed, comparing the 
measured mass distributions in all reconstructed channels to their simulated counterparts. 

Some small differences are observed in the 3h and /i27r° channels. However no 
new evaluation of this systematic uncertainty was done for this analysis as it is not 
dominant. The result is taken from the published paper on 1991-1993 data [3] which 
was based on detailed tests of several models describing the Oi resonance. Similarly, 
systematic uncertainties were derived for the higher multiplicity modes by comparing 
decay models with different resonance contributions in the final states. Finally the effect 
of T polarization on the efficiency matrix is estimated to be negligible. 

The corresponding errors are listed in Table 11 in the column 'dyn'. 



10.16 Monte Carlo statistics 

Multinomial fluctuations for each generated class in the Monte Carlo sample using the 
measured fractions of events in each reconstructed class are included in the efficiency 
matrix uncertainties. Generating 1000 sets of efficiency matrices, the RMS of the 
branching ratio values are obtained and taken as the uncertainty from the finite Monte 
Carlo statistics. It is shown in Table 11 in the column 'mcs'. 
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10.17 Total systematic errors 



Assuming that all the above systematics from specific sources are uncorrelated, the total 
systematic error is computed for each channel. The final values are listed in the last 
columns of Table 11. Correlations exist between different channels for a given source 
of systematics: the correlation matrix for the total systematic uncertainties is given in 
Section 13.3. 



11 Global systematic checks 

After taking into account all the corrections to the simulation (measured r branching ratio 
values, fake photon simulation, . . . ) and inchiding non-r backgrounds, some distributions 
are compared in data and Monte Carlo in order to provide global checks of the consistency 
between data and simulation. Of course such tests cannot be taken hterally as they could 
reveal some difference in the r decay dynamics in some channel, however such effects 
are expected to be small as most r decay channels are well understood from the basic 
principles of the Standard Model. This comparison is done for the distributions of charged 
particle momenta, of photon and tt^ energies and of the hadronic invariant mass spectra 
for each r decay mode. 

11.1 Charged particle momentum distributions 

The only accessible observable in the leptonic channels is the lepton momentum. 
Figures 31 and 32 show the respective track momentum distributions for electrons and 
muons. They are in good agreement with the simulation including backgrounds and the 
Standard Model spectrum. 

The charged hadron momentum distribution is shown in Fig. 33. Good agreement is 
observed between data and the simulation. 

11.2 Photon energy distribution 

The comparison of photon energy spectra in data and Monte Carlo is shown in Figs. 34 
and 35 for 1991-1993 and 1994-1995 data sets respectively. Good agreement is observed 
except for the largest energies. This effect is understood as coming from the correction of 
saturation and leakage which is applied to single clusters in data, while no such correction 
is done in the simulation. This procedure is incorrect for unresolved 7r°'s not recognized 
as such, since the correction is applied to the resulting shower rather than to the two 
merged individual photonic showers, leading to an overestimate of the 7r° energy. This 
shortcoming has no effect on the branching ratio measurement. The agreement between 
data and Monte Carlo at low energy after fake photon correction shows that the procedure 
is reliable, since no energy information was used in the fake photon correction procedure. 
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Figure 31: Comparison of electron momentum spectra in data (points) and simulation 
(histogram) for the full 1991-1995 data sample. The shaded histogram is the contribution 
of non-r background. 
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Figure 32: Comparison of muon momentum spectra in data (points) and simulation 
(histogram) for the full 1991-1995 data sample. The shaded histogram is the contribution 
of non-r background. 
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Figure 33: Comparison of charged hadron momentum spectra in data (points) and 
simulation (histogram) for the full 1991-1995 data sample. The shaded histogram is 
the contribution of non-r background. 
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Figure 34: Comparison of photon energy spectra for all r events in data (points) and 
simulation (histogram) after all the corrections for 1991-1993 data set. The dotted 
line shows the simulated distribution before correction and the shaded histogram is the 
contribution of non-r background. 
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Figure 35: Comparison of photon energy spectra for all r events in data (points) and 
simulation (histogram) after all the corrections for 1994-1995 data set. The dotted 
line shows the simulated distribution before correction and the shaded histogram is the 
contribution of non-r background. 
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11.3 TT energy distribution 



The comparison of vr" energy spectra in data and Monte Carlo is shown in Figs. 36 and 37 
for 1991-1993 and 1994-1995 data sets respectively. Good agreement is observed except 
at energies below 1 GeV for the 1994-1995 data and for the highest energies. The former 
discrepancy is within the quoted systematic uncertainty for 7 — tt" reconstruction. As 
for high energy, the explanation is the same as for single photons since, at high energy, 
unresolved 7r°'s dominate the sample. 

11.4 Invariant mass spectra in multihadron channels 

The comparison of mass spectra for all r decays in data and Monte Carlo is shown in 
Fig. 38 for the full 1991-1995 data set. As noticed earlier perfect agreement cannot be 
expected since the Monte Carlo simulation of the dynamics in some decay modes may not 
be accurate. Even in prominent channels like tttt^, Sti and 7r27r°, the precise shape of the 
mass distribution is not predicted by theory and the phenomenology of the corresponding 
spectral functions docs, on the contrary, rely on precise measurements (see Section 15 on 
the measurement of spectral functions from this analysis). The long tail above the r mass 
in the distribution is dominated by secondary interaction of hadronic tracks with inner 
detector material. Globally, the agreement is reasonable and confirms that the simulation 
of secondary interactions in the detector is adequate. 

The comparison of mass spectra for reconstructed hn'^ events in data and Monte 
Carlo is shown in Fig. 39 for the full 1991-1995 data set, with the Monte Carlo predicted 
feedthrough from other channels and non-r backgrounds. 

12 From reconstructed classes to exclusive modes 

So far branching fractions have been determined in 13 classes corresponding to major r 
decay modes. However, as shown in Table 3, these classes still contain the contributions 
from final states involving kaons. The latter are coming from Cabibbo-suppressed r decays 
or modes with a KK pair, both characterized by small branching ratios compared to the 
nonstrange modes without kaons. 

Complete analyses of r decays involving neutral or charged kaons have been performed 
by ALEPH on the full LEP 1 data [6, 7, 8]. They are summarized in Ref. [9] where 
measurements with Kg or are combined. The ALEPH analyses have provided 
measurements of branching ratios of modes with kaons containing up to 4 hadrons in 
the final states. Thus they are fully adequate to cover the needs of the present analysis 
of the nonstrange modes. 

The r decays involving r] or u mesons also require special attention in this analysis 
because of their electromagnetic decay modes. Indeed the final state classification relies in 
part on the 7r° multiplicity, thereby assuming that all photons — except those specifically 
identified as bremsstrahlung or radiative — originate from tt^ decays. Therefore the non- 
7r° photons from rj and cu decays are treated as tt^ candidates in the analysis and the 
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Figure 36: Comparison of 7r° energy spectra for all r events in data (points) and simulation 
(histogram) after all the corrections for 1991-1993 data set. The dotted line shows the 
simulated distribution before correction and the shaded histogram is the contribution of 
non-r background. 
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Figure 37: Comparison of vr" energy spectra for all r events in data (points) and simulation 
(histogram) after all the corrections for 1994-1995 data set. The dotted line shows the 
simulated distribution before correction and the shaded histogram is the contribution of 
non-r background. 
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Figure 38: Comparison of mass spectra (assuming all charged particles are pions) in data 
(points) and simulation (histogram) after fake photon correction for all the hadronic r 
decays (except single hadron) in 1991-1995 data sample. The shaded histogram is the 
contribution of non-r background. The same plots are displayed in logarithmic (top) and 
linear (bottom) vertical scales. 
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Figure 39: Comparison of hadronic mass spectra in data (points) and simulation 
(histogram) after fake photon correction for the /ivr^sample in 1991-1995 data. The lower 
shaded histogram is the contribution of non-r background and the upper shaded is from 
T feedthrough. The plots are displayed in logarithmic (top) and linear (bottom) vertical 
scales. 
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Table 12: Corrections for the exclusive nonstrange branching ratios. 'QE' and 'E' denote 
the quasi-exclusive (with kaons, uo and r] included) and exclusive modes, respectively. 
Treating the uj and r} contributions separately is made necessary because of their large 
radiative {i.e. with photons not originating from 7r"'s) modes. 



QE class E class correction to BR (%) 



e 


e 


-0.000 


± 


0. 


,000 






-0.000 


± 


0. 


,000 


h 


TT" 


-1.341 


± 


0. 


,040 


h7T° 




-0.756 


± 


0. 


,038 


h2TT^ 


7r-27r° 


-0.408 


± 


0. 


,030 






-0.236 


± 


0, 


,032 


/i47rO 




-0.085 


± 


0. 


,016 


3h 


7r~7r~7r''" 


-0.770 


± 


0. 


,057 


3/i7r° 


7r~7r~7r'^7r° 


-1.994 


± 


0. 


,100 




7r~7r~7r"'"27r° 


-0.480 


± 


0. 


,071 




7r~7r~7r"'"37r° 


-0.032 


± 


0. 


,006 


5h 


37r-27r+ 


-0.026 


± 


0. 


,004 


5/i7r° 


37r-27r+7r° 


-0.012 


± 


0. 


,002 



systematic bias introduced by this effect must be evaluated. The corrections are based 
on specific measurements by ALEPH of r decay modes containing those mesons [13]. 
Thus the final results correspond to exclusive branching ratios obtained from the values 
measured in the topological classification, corrected by the removed contributions from 
K, 1] and uj modes measured separately, taking into account through the Monte Carlo 
their specific selection and reconstruction efficiencies to enter the classification. This 
delicate bookkeeping takes into account all the major decay modes of the considered 
mesons [25], including the isospin- violating cj — > 7r"'"7r~ decay mode. The main decay 
modes considered arc nuj, nn'^uj and nn'^ri with branching fractions of (2.26±0.18) x 10~^, 
(4.3 ± 0.5) X 10"■^ and (1.80 ± 0.45) x 10"^ [13], respectively. The first two values are 
derived from the branching ratios for the 37r7r'' and 37r27r° modes obtained in this analysis 
and the measured u; fractions of 0.431 ± 0.033 from ALEPH [13] and the average value, 
0.78 ± 0.06, from ALEPH [13] and CLEO [26], respectively. 

Some much smaller contributions with 77 and cu have been identified and measured by 
CLEO [27] with the decay modes r lyrrjn-n+TT- ((2.4 ± 0.5) x lO""^), r 1/^7777-277° 
((1.5 ± 0.5) X 10"^), r u^LUTT-TT+TT- ((1.2 ± 0.2) X 10-^), and r ^ z/^o^7r-27r° 
((1.5 ±0.5) X lO^^). Even though the corrections from these channels are very small they 
have been included for the sake of completeness. Finally, another very small correction 
has been applied to take into account the ai radiative decay into ttj with a branching 
fraction of (2.1 ± 0.8) x 10"^ obtained from Ref. [28]. 

The corrections used to obtain exclusive branching ratios for the listed nonstrange 
modes are given in Table 12. 
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13 Results 



13.1 Overall consistency test 

Rejected r hemispheres because of charged particle identification cuts are placed in class 
14; these cuts include the 2 GeV minimum momentum and the ECAL-crack veto for some 
one-prong modes, and the strict definition of higher multiplicity channels. As already 
emphasized, this sample does not correspond to a nominal r decay mode and should be 
explained by all other measured fractions in the other classes and the efficiency matrix. 
Thus the determination of a hypothetical signal in this class is a measure of the level of 
consistency achieved in the analysis. 

For this determination the efficiency of the possible signal in class 14 is taken to be 
100%. The results, already shown in Table 8 separately for the 1991-1993 and 1994- 
1995 data sets, are consistent and are combined to give B14 = (0.066 ± {0.027)stat =t 
{0.021)syst,c i (0.025)syst,u„c)%, where the last two errors refer to the common and 
uncommon uncertainties from the two data sets. With a combined error of 0.042% this 
value is consistent with zero and provides a nontrivial check of the overall procedure at 
the 0.1% level for branching ratios. It is interesting to note that this value coincides, 
approximately and accidentally, with the limit achieved of 0.11% at 95% CL in a direct 
search for "invisible" decays not selected in the 13-channel classification. 

In the following it is assumed that all r decay modes have been properly considered at 
the 0.1% precision level and no physics contribution beyond standard r decays is further 
allowed. Thus the quantity B14 is now constrained to be zero. 

It can be further noticed that this analysis provides a branching ratio in the SttStt^ class 
which is consistent with zero for both 1991-1993 and 1994-1995 data sets (see Table 8). 
The result is therefore given as an upper limit at 95% CL 

^3,3.0 < 4.9 X 10-^ (15) 

consistent with the measurement made by CLEO [29] yielding B^t^^t^o = (2.2 ±0.5) x 10~^. 
The final state is dominated by rj and u resonances [29] and using other channels allows 
a lower limit to be obtained for this branching ratio, (2.6 ± 0.4) x 10~^. In the following 
a value of (3 ± 1) x 10~^ is used as input for this channel and the global analysis is 
performed in terms of the remaining 12 defined channels which are refitted. As for other 
channels proper subtractions are made for the contributions of modes with 77 which are 
listed separately. 

13.2 Comparison of 1991-1993 and 1994-1995 results 

Since the same procedure is applied for the analyses of 1991-1993 and 1994-1995 data the 
results must be consistent within the statistical errors of data and Monte Carlo. Table 13 
shows the fist of the differences of branching ratios with their expected fluctuations. Good 
agreement is observed with a of 8.6 for 11 DP. 



70 



Table 13: Differences of branching ratios between 1991-1993 and 1994-1995 data samples; 
only the statistical errors in data and Monte Carlo are considered. 



class 


A BR (%) 


e 


0.040 ± 0.148 




0.061 ± 0.140 


h 


0.145 ± 0.139 


7 n 


-0.186 ± 0.197 


/i27r° 


0.018 ± 0.178 


/i37r° 


-0.105 ± 0.143 


M7r° 


0.065 ± 0.076 


3h 


0.163 ± 0.128 


3/i7rO 


-0.201 ± 0.120 


3/i27rO 


0.046 ± 0.078 


5h 


0.006 ± 0.018 


5hn° 


-0.006 ± 0.014 



In conclusion the two independent data and Monte Carlo samples give consistent 
results. The 1994-1995 results confirm the trend of larger hir^ and Shn^, and smaller h 
and 3h branching ratios compared to the previous analyses, as observed in 1991-1993 data 
sample. 

13.3 Final combined results 

Finally the two sets of results are combined. Using only statistical or total weights — in the 
latter taking into account correlated errors from dynamics and secondary interactions — 
gives almost identical results. The final results obtained with the total weights are shown 
in Table 14. 

The branching ratios obtained for the different channels are correlated with each 
other. On one hand the statistical fiuctuations in the data and the Monte Carlo sample 
are driven by the multinomial distribution of the corresponding events, producing well- 
understood correlations. On the other hand the systematic effects also induce significant 
correlations between the different channels. All the systematic studies were done keeping 
track of the correlated variations in the final branching ratio results, thus allowing a 
proper propagation of errors. The full covariance matrices from statistical (from data) 
and systematic origins are given in Tables 15 and 16. 

13.4 An independent analysis of the leptonic branching ratios 

A dedicated analysis of the leptonic branching ratios has been performed as a by-product 
of one of the two methods used by ALEPH to measure the r polarization and polarization 
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Table 14: Combined results for the exclusive branching ratios (B) for modes without 
kaons. The contributions from channels with 77 and uj are given separately, the latter 
only for the electromagnetic cu decays. All results are from this analysis, unless explicitly 
stated. The "estimates" are discussed in the text. 



mode 
e 

TT" 

7r-27r° 
7r-37r° 
7r-47r° 
7r~7r~7r"' 

7r~7r~7r"'"27r° 
7r~7r~7r'^37r° 
37r"27r+ 
37r-27r+7r° 

7r~27r°77 

ar(- 



B it^stat ± 



syst 



TT 7) 



7r°7, TT+TT" 



7r~7r°a;(— > 7r°7, tt+tt^) 
7r~27r°a;(— >• 7r°7, 7r"'"7r~) 
7r~7r~7r''"a;(— > 7r°7, 7r''"7r~) 



17 


837 


-1- 


n 


072 




n 

yj 


036 




17 


319 


± 





070 


± 





032 




10 


828 


± 





070 


± 





078 




25 


471 


± 





097 


± 





085 




9 


239 


± 





086 


± 





090 







977 







069 


± 





058 







112 


± 





037 


± 





035 




9 


041 


± 





060 


± 





076 




4 


590 


± 





057 


± 





064 







392 


± 





030 


± 





035 







013 


± 





000 


± 





010 


estimate 





072 


± 





009 


± 





012 







014 


± 





007 


± 





006 







180 


± 





040 


± 





020 


ALEPH [13] 





015 


± 





004 


± 





003 


CLEO [27] 





024 


± 





003 


± 





004 


CLEO [27] 





040 


± 





000 


± 





020 


estimate 





253 


± 





005 


± 





017 


ALEPH [13] 





048 


± 





006 


± 





007 


ALEPH [13] + CLEO 





002 


± 





001 


± 





001 


CLEO [27] 





001 


± 





001 


± 





001 


CLEO 27 
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Table 15: Correlation matrix of the statistical errors on the branching fractions. 







h 




/i27r" 






■ih 


3hTT^ 


3/i27r" 


3/i37r" 


5h 


5/i7r" 


e 


-0.21 


-0.15 


-0.25 


-0.09 


-0.01 


0.00 


-0.15 


-0.10 


0.03 


-0.06 


0.00 


0.01 




1.00 


-0.13 


-0.21 


-0.07 


-0.06 


0.00 


-0.09 


-0.07 


0.00 


-0.02 


0.00 


-0.04 


h 




1.00 


-0.31 


-0.02 


0.01 


-0.06 


-0.12 


-0.06 


-0.02 


0.01 


-0.01 


0.02 


hTT^ 






1.00 


-0.40 


0.05 


0.00 


-0.11 


-0.06 


-0.02 


0.00 


-0.04 


-0.04 










1.00 


-0.51 


0.26 


-0.09 


0.01 


-0.07 


0.06 


-0.01 


0.03 












1.00 


-0.75 


0.01 


-0.03 


0.05 


-0.02 


-0.01 


0.01 


Mtt^ 












1.00 


-0.02 


-0.02 


-0.03 


0.01 


0.02 


-0.03 


3h 














1.00 


-0.33 


0.08 


-0.05 


-0.04 


0.00 


3/i7r° 
















1.00 


-0.45 


0.19 


-0.02 


-0.02 


3/i27r° 


















1.00 


-0.65 


0.03 


0.02 


3/i37r° 




















1.00 


-0.01 


-0.04 


5h 






















1.00 


-0.24 


























1.00 



Table 16: Correlation matrix of the systematic errors on the branching fractions. 







h 




/i27r" 


/i37r" 


/i47r" 


3/i 


3/i7r" 


3/i27r" 


3/i37r" 


hh 


5/i7r" 


e 


-0.17 


-0.01 


0.02 


0.01 


0.03 


-0.08 


-0.17 


-0.22 


-0.05 


0.02 


0.00 


0.00 




1.00 


0.05 


0.09 


-0.03 


0.02 


-0.13 


-0.11 


-0.24 


-0.06 


0.01 


0.03 


-0.04 


h 




1.00 


0.36 


-0.29 


-0.32 


-0.42 


0.34 


-0.40 


-0.40 


-0.07 


0.16 


-0.09 








1.00 


-0.35 


-0.02 


-0.33 


0.01 


-0.54 


-0.26 


0.02 


0.11 


-0.06 


/i27rO 








1.00 


-0.01 


0.13 


-0.24 


0.07 


0.13 


0.06 


-0.13 


0.03 


/i3^o 










1.00 


-0.13 


-0.29 


-0.02 


0.15 


0.09 


-0.06 


0.04 


/i47r° 












1.00 


-0.14 


0.34 


0.27 


0.00 


-0.12 


-0.05 


3h 














1.00 


-0.03 


-0.16 


-0.11 


0.17 


-0.06 


3/i7r° 
















1.00 


0.04 


-0.03 


-0.07 


-0.01 


3/i27r0 


















1.00 


-0.14 


-0.09 


0.07 


3/i37r° 




















1.00 


-0.02 


0.02 


5h 






















1.00 


-0.26 


























1.00 
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asymmetry [16]. It identifies and selects directly single r hemispheres with leptons without 
a full selection of other tau decay modes. In order to normalize the selected samples, the 
number of produced r pairs is derived using the e'^e" — > t"^t~ cross section measured 
by ALEPH [19] and the precise determination of luminosity using small-angle Bhabha 
scattering. 

Particle identification is performed using a method quite similar to the one used in the 
global analysis and described in Section 4.2 but using a totally independent code. The 
performances of the two procedures are very close. Similarly, the treatment of photons 
and the separation between genuine and fake photons follows similar philosophies, but it 
is implemented separately. 

Hemispheres with r decay leptonic candidates are selected in two steps. In the 
preselection step, events with acollinear "jets" (cos^acoi < —0.9) are retained with a loose 
lepton identification at least on one side. Then, strict cuts are applied to select electron 
and muon hemispheres, by rejecting hadronic r decays (electron or muon identification, 
veto if a 7r° is reconstructed in the same hemisphere) and non-r background (using 
methods similar to those described in Section 3.3). The selection efficiencies are about 
69% and 75% for electronic and muonic decays, respectively. The corresponding values 
for the r decay fccdthrough contaminations are 1.0% and 0.8%, and 1.5% and 0.4% for 
the non-r background, respectively. The muon identification in this method allows the 
detection of lower muon momenta down to 1.3 GeV, compared to 2 GeV in the global 
method. A total of 48882 electron and 50782 muon hemispheres are thus selected. In this 
method the largest systematic uncertainty originates from the r-pair normalization and 
is completely correlated for the electron and muon channels. The results of this specific 
analysis are: 

B{T^VreV^) = (17.778 ± 0.080 ± 0.049)% , (16) 
B{t iyrfJ''J^n) = (17.299 ± 0.077 ± 0.045)% , (17) 

where the first errors are statistical and the second systematic. 

Although both analyses share the same data and have both large selection efficiencies 
they have been carried out in a completely different and independent way, from the 
selection of candidates to the identification method of the final state. The comparison 

thus provides a valuable check. 

The samples of events have been compared and the correlations found to be 0.861 for 
electrons and 0.896 for muons. Part of the Monte Carlo statistics is also non common. 
An important effect could come from the normalization of the independent method using 
the luminosity-derived number of produced rr events, which is common to both electron 
and muon analyses. 

The results of the independent analysis are consistent with the present values of the 
leptonic branching ratios. The differences found, AB^ = (0.059 ± 0.051 ± 0.029))% 
and AB^ = (0.020 ± 0.042 ± 0.028)%, are well within the expectation of the respective 
noncommon errors. The second quoted errors are from the normalization of rr events 
and are completely correlated. 
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Since the results are in agreement they could be combined. However, the resulting 
improvement is not visible for the statistical error and small for the (non-dominant) 
systematic uncertainty. Furthermore, it is very cumbersome to average the leptonic 
branching fractions alone in the framework of the global analysis where all branching 
ratios are derived in a consistent way. For these reasons the final results are taken from 
the global analysis, while the independent analysis provides a meaningful check for the 
leptonic channels. 



14 Discussion of the results 

14.1 Comparison with other experiments 

As the results from other experiments are most often given without K — tt separation 
for charged hadrons, comparison is made summing branching ratios with charged kaons 
and pions. Also contributions from modes with rj are included in the relevant hadronic 
final states. Figures 40, 41, 42, 43, 44, 45, 46, 47, 48, 49 and 50 show that the results 
of this analysis are in good agreement with those from other experiments. In all these 
cases, ALEPH achieves the best precision, except for branching ratios under 0.1%, such 
as r — > 5hu and r — > Shn^u, where CLEO results are more precise because of their higher 
statistics. 

A meaningful comparison can be performed between the exclusive fractions and the 
topological branching ratios Bi, where i refers to the charged particle multiplicity in the 
decay. Even though the latter have essentially no physics interest, their determination 
can constitute a valuable cross check as they depend only on selection efficiency, tracking, 
handling of secondary interactions and electron identification for photon conversions, and 
not on photon identification. The results from this analysis can be compared in this way 
with a dedicated analysis recently performed by DELPHI [48]. Summing up appropriately 
(both analyses assume a negligible contribution from hadronic multiplicities larger than 
five, in agreement with the 90% CL hmit by CLEO [49], Bj < 2.4 x 10"^), one gets 



in good agreement with the DELPHI values, B3 = (14.569 ± 0.093 ± 0.048)% and 
B5 = (0.115 ± 0.013 ± 0.006)%. The rather small systematic uncertainty in the DELPHI 
results reflects the fact that a sharper study of hadronic interactions can be performed 
when only charged particles are considered in the analysis. In addition the modes with 
Kg — > TT'^TT" decays are subtracted statistically here, rather than trying to identify them 
on an event-by-event basis. 



S3 = (14.652 ± 0.067 ± 0.086)% , 
= (0.093 ± 0.009 ± 0.012)% , 



(18) 
(19) 
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Figure 40: Comparison of ALEPH measurement with published precise resuhs from other 
experiments for r evv. References for other experiments are ARGUS [30], CLEO [31], 
DELPHI [32], OPAL [33], L3 [34]. 
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Figure 41: Comparison of ALEPH measurement with published precise results from other 
experiments for r — > ^vv. References for other experiments are ARGUS [30], CLEO [31], 
DELPHI [32], OPAL [35], L3 [34]. 
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Figure 42: Comparison of ALEPH measurement with published precise resuhs from other 
experiments for r ^ hv (sum of and Kv). References for other experiments are 
CLEO [31], OPAL [36]. 
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Figure 43: Comparison of ALEPH measurement with published precise results from other 
experiments for r h-K^v (sum of titx^v and Ktt^u). References for other experiments 
are CLEO [37], OPAL [36]. 
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Figure 44: Comparison of ALEPH measurement with published precise resuhs from other 
experiments for r h2Tx^v (sum of it2tt^i' and K2t{^v). References for other experiments 
are CLEO [38]. 
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Figure 45: Comparison of ALEPH measurement with published precise results from other 
experiments for r — > /i(3,4)7r°z/ (sum of 7r(3,4)7r°z/ and 7^(3, 4)7r°z/). References for other 
experiments are CLEO [38]. 
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Figure 46: Comparison of ALEPH measurement with published precise results from 
other experiments for r Shu (sum of Svrz/, K2tiv and 2KTiv). References for other 
experiments are CLEO [39], OPAL [40]. 
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Figure 47: Comparison of ALEPH measurement with published precise results from other 
experiments for r Sh-K^u (sum of Stttt^u, K2tiv and 2Kttv). References for other 
experiments are CLEO [39]. 
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Figure 48: Comparison of ALEPH measurement with published precise resuhs from other 
experiments for r 3/i(2,3)7rV (sum of 37r(2,3)7rV, K2tt{2,3)tt^u and 2Js:7r(2, 3)7rV). 
References for other experiments are CLEO [41]. 
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Figure 49: Comparison of ALEPH measurement with published precise results from other 
experiments for r dhu. References for other experiments are HRS [42], ARGUS [43], 
CLEO [44], OPAL [45], BABAR [46]. 
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Figure 50: Comparison of ALEPH measurement with published precise results from other 
experiments for r 5/i7r°z/. References for other experiments are CLEO [47], OPAL [45]. 
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14.2 Universality in the leptonic charged current 



14.2.1 ji — e universality from the leptonic branching ratios 

In the standard V-A theory with leptonic couphng gi at the WlVi vertex, the r leptonic 
partial width can be computed, including radiative corrections [50] and safely neglecting 
neutrino masses: 

r(.^.,i.,(T))^^/(^)^^^;, (20) 



where 

Gi = 



9l 



5^ - 1 + ^^' 



7 



f{x) = 1 - 8x + - - 12x^lnx . (21) 

Numerically, the W propagator correction and the radiative corrections are small: 5^ = 
1 + 2.9 X 10-^ and (5; = 1 - 43.2 x 10"! 

Taking the ratio of the two leptonic branching fractions, a direct test oi /i — e 
universality is obtained. The measured ratio 

= 0.9709 ± 0.0060 ± 0.0029 (22) 

Be 

agrees with the expectation of 0.97257 when universality holds. Alternatively the 
measurements yield the ratio of couplings 

^ = 0.9991 ± 0.0033 (23) 

ge 

which is consistent with unity. 

This result is in agreement with the best test ol ji — e universality of the W couplings 
obtained in the comparison of the rates for tt — ^ /xI7^ and tt — > eVf. decays, where 
the results of the two most accurate experiments [51, 52] can be averaged to yield 
— = 1.0012 ± 0.0016. The results have comparable precision, but it should be pointed 
out that they are in fact complementary. The r result given here probes the coupling to a 
transverse W (hehcity ±1) while the tt decays measure the coupling to a longitudinal W 
(hclicity 0). It is indeed conceivable that either approach could be sensitive to different 
nonstandard corrections to universality. 

Since and are consistent with fi — e universality their values can be combined, 
taking common errors into account, into a consistent leptonic branching ratio for a 

2 

massless lepton (the electron, since /(^) = 1 within 10~^) 

Qimi=o) ^ ± 0.044 ± 0.022)% , (24) 

where the first error is statistical and the second systematic. 
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14.2.2 Tests of r — and r — e universality 

Comparing the rates for r(r — > VT-eV , V{t — > I'rlM'ij.il)) and r(// — > f^ei7e(7)) 
provides direct tests of the universahty of r — fi — e couplings. Taking the relevant 
ratios with calculated radiative corrections, one obtains 

f^V = ^f^)^e^A^A,, (25) 



9t 



2 



5^— -f-A,^A^ , (26) 



where /(^) = 0.9998, A^^ = ^ = 1 - 2.9 x 10-^ A^ = ^ = 1 + 8.5 x 10-^ and is 
the lepton I lifetime. 



From the present measurements of i?e, -B^, the r mass [25], = (1777.03^^:^1]) MeV 
(dominated by the BES result [66]), the r lifetime [25], = (290. 6± 1.1) fs and the other 
quantities from Ref. [25], universality can be tested: 

— = 1.0009 ± 0.0023(Se) ± 0.0019(t^) ± 0.0004(m^) , (27) 

— = 1.0001 ± 0.0022(5J ± 0.0019(r^) ± 0.0004(m^) . (28) 

9e 



14.2.3 T — fjL universality from the pionic branching ratio 

The measurement of also permits an independent test of r — universality through 
the relation 



9r\ T^2m^m^ / I 



1 - 4 i ^r/^ 



(29) 



where the radiative correction [53] amounts to 5,-/^ = 1.0016 ± 0.0014. Using the world- 
averaged values for the r and tt (r^) lifetimes, and the branching ratio for the decay 
TT —>■ iiu [25], the present result for B^^, one obtains 

— = 0.9962 ± 0.0048(5^) ± 0.0019(r^) ± 0.0002(rrv) ± 0.0007(rarf) , (30) 
9fi 

comparing the measured value, = (10.828 ± 0.105)%, to the expected one assuming 
universality (10.910 ± 0.042)%. 

The two determinations of ^ obtained from Be and B^^ are consistent with each other 

9iJ. 

and can be combined to yield 

— = 1.0000 ± 0.0021(5e, B^) ± 0.0019(r^) ± 0.0004(m0 . (31) 
9iJ- 
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Figure 51: The measured value for Be compared to predictions from other measurements 
assuming leptonic universahty. The vertical band gives the average of all determinations. 

Universality of the r and fi charged-current couplings holds at the 0.29% level with about 
equal contributions from the present determination of Bg and B^, and the world- averaged 
value for the r lifetime. 

The consistency of the present branching ratio measurements with leptonic universality 
is displayed in Fig. 51 where the result for B^ is compared to computed values of B^ using 
as input B^ (assuming e — /i universality), r,- and (/U — t universality), and and 
{fi — T universality). All values are consistent and yield the average 

^universality _ _QIQ ± 0.039)% . (32) 



14.3 The TTTT branching ratio in the context of 

The tttt" final state is dominated by the p resonance as demonstrated in Fig. 39. Its mass 
distribution — or better, the corresponding spectral function, see Section 15 — is a basic 
ingredient of vacuum polarization calculations, such as used for computing the hadronic 
contribution to the anomalous magnetic moment of the muon a^"''^. In this case the p 
contribution is dominant (71%) and therefore controls the final precision of the result. It 
was proposed in Ref. [54] to use the spectral functions obtained from the measurement 
of hadronic r decays in order to improve the precision of the prediction for a'^f'^. The 
calculation was later improved with the help of QCD constraints for energies above the r 
mass [55] and even below [56]. 

The normalization of the spectral function is provided by the branching fraction -Btt^o . 
The present world average is completely dominated by the published ALEPH result [3]. 



89 



nn 



-0.6 -0.4 -0.2 
B 



0.2 0.4 0.6 



measured predicted 



(%) 



Figure 52: The differences between the measured values for B.,^ and Bt^2tt0 and their 
respective predictions from leptonic universahty and isospin symmetry with vr-mass 
breaking, compared to the precision on Bt,t,o. 



The new result given here is larger by 0.68%, thus one can expect a slightly larger 
contribution to a!^'^'^. 



A new evaluation [57] of a|^'^'^ was available, using the preliminary spectral functions 
from the present analysis, the published CLEO results [58] and new results from e^e~ 
annihilation from CMD-2 [59]. Revision of the CMD-2 results [60] prompted a re- 
evaluation [61], which revealed a disagreement between the r and e~^e~ spectral functions 
(see Section 15). Whereas the r estimate leads to a prediction of the muon magnetic 
moment in agreement with the latest most precise measurement from the BNL experiment 
E-821 [62], the predicted value using only e^e~ data lies 2.4 standard deviations below 
the measurement. In view of this situation, it is important to check all the ingredients, 
in particular the determination of the branching ratio Bj^^^o. As most of the systematic 
uncertainty in Bt^^o comes from 7/71° reconstruction, it is helpful to cross check the results 
in the "adjacent" hadronic modes, i.e. the vr and 7r27r° channels. This is possible if 
universality in the weak charged current is assumed, leading to an absolute prediction of 
Bt, using as input the r and n lifetimes (Section 14.2.3, Eq. (29)), and by computing Bt,2tt0 
assuming isospin symmetry from the measurement of B^^r which is essentially uncorrelated 
with -BttttO (Section 14.4). The comparison, shown in Fig. 52, does not point to a systematic 
bias in the determination of B^^^o within the quoted uncertainty. 
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14.4 ai decays to Sir and 7r27r' 



With the level of precision reached it is interesting to compare the rates in the Stt and 

7r27r° channels which arc dominated by the Oi resonance. The dominant pn intermediate 
state leads to equal rates, but a small isospin-breaking effect is expected from different 
charged and neutral vr masses, slightly favouring the 7r27r° channel, as discussed below. 

A recent CLEO partial-wave analysis of the 7r27r'' final state [63] has shown that the 
situation is in fact much more complicated with many intermediate states, in particular 
involving isoscalars, amounting to about 20% of the total rate and producing strong 
interference effects. A good description of the Oi decays was achieved in the CLEO study, 
which can be applied to the Sir final state, predicting [63] a ratio of the rates 3'it/tt2tt^ 
equal to 0.985. This value, which includes known isospin-breaking from the pion masses, 
turns out to be in good agreement with the measured value from this analysis which shows 
the expected trend 

— ^ = 0.979 ± 0.018 . (33) 

-D7r27rO 



14.5 Separation of vector and axial-vector contributions 
14.5.1 Strange and nonstrange contributions 

From the complete analysis of the r branching ratios presented in this paper, it is possible 
to determine the nonstrange vector (V) and axial-vector (A) contributions to the total 
r hadronic width, conveniently expressed in terms of their ratios to the leptonic width, 
called Rry and Rr,A, respectively. The determination of the strange counterpart Rt,s is 
already pubhshed [9] . 

The ratio Rr for the total hadronic width is calculated from the leptonic branching 
ratios alone, and eventually from the electronic branching ratio alone, 

Rr = ^ ^ = — - 1.9726 

Be -Be 

= 3.642 ± 0.012 . (34) 

taking for B{t~ — > e~ Pe^r) the value obtained in Section 14.2 assuming universality in 
the leptonic weak current. All R values given below are rescaled so that the sum of the 
hadronic and the leptonic branching ratios, the latter computed using the "universal" 
value given in Eq. 32), add up to unity. The deviations introduced in this way are 
extremely small, less than one tenth of the experimental error, but this procedure 
guarantees the consistency between all values. Using the ALEPH measurement of the 
strange branching ratios [9], supplemented by the small contribution from the K*~r) 
channel measured by CLEO [64], 

Bs = (2.85 ± 0.11) % , 
Rr,s = 0.1603 ± 0.0064 , (35) 
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the following result is obtained for the nonstrange component 

Bv+A = (62.01 ± 0.14) % , 
R^y+A = 3.482 ± 0.014 . (36) 

Separation of V and A components in hadronic final states with only pions is 
straightforward. Isospin invariance relates the spin-parity of such systems to their number 
of pions: G-parity =1 (even number) corresponds to vector states, while G = —1 (odd 
number) tags axial-vector states. This property places a strong requirement on the 
efficiency of 7r° reconstruction, a constraint that was emphasized in this analysis. 

14.5.2 V and A separation in modes with a KK pair 

Modes with a KK pair are not in general eigenstates of G-parity and contribute to both 

V and A channels. While the decay to K~K^ is pure vector, the situation is a priori not 
clear in the KKn channel, observed in three charged modes: K^K^tt~, K~K^ti^ and 
K^K 7r~. Three sources of information exist on the possible V / A content in this decay 
mode: 

(i) in the ALEPH analysis of r decay modes with kaons [9], an estimate of the vector 
contribution was obtained using the available e+e~ annihilation data in the KKtt channel, 
extracted in the 1=1 state. In fact, this contribution was found to be very small, yielding 
a branching ratio consistent with zero, Bcyc{t Ur[KKTi]v) = (0.26 ± 0.39) x 10~^, 
corresponding to an axial fraction, /^(CVC) = 0.941o;q|. 

(ii) The CVC/ ALEPH result is corroborated by a partial- wave and lineshape analysis 
of the ai resonance from r decays in the ^'i-7r~27r'^ mode by CLEO [63]. The observation 
through unitarity of the opening of the K*K decay mode of the ai is claimed and a 
branching ratio B{ai K*K) = (3.3 ± 0.5)% obtained. Knowing the decay rate of 
T —>■ i'tO'I, it is easy to sec that such a result completely saturates, and even exceeds, 
the total rate for the KKn channel. The corresponding axial fraction turns out to be 
/a(CLEO 37r) = 1.30 ±0.24. 

(in) A new piece of information, also contributed by CLEO, but conflicting with 
the two previous results, was recently published [65]. It is based on a partial- wave 
analysis in the K^K^ir^ channel using two-body resonance production and including 
many possibly contributing channels. A much smaller axial contribution is obtained, 
/a(CLEO KKn) = 0.56 ± 0.10. 

The three determinations are displayed in Fig. 53. Since they are not consistent, a 
conservative value of Ja = 0.75 ± 0.25 is assumed. It encompasses the range allowed by 
the previous results and still represents some progress over our former analyses [11, 12] 
where a value of 0.5 ±0.5 was used. For the decays into KKnn no information is available 
in this respect and the same conservative fraction of 0.5 ± 0.5 is assumed. 

All the present measurements and previous ones for modes with a KK pair are 
summarized in a complete list of determined branching ratios for nonstrange hadronic 
final states up to six hadrons, and four with kaons. The results are given in Table 17. 
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Figure 53: The axial- vector fraction J'a of tlie liadrons in tlie r — > UrKKir decay mode 
from three independent determinations: ALEPH [9] using CVC and e~^e~ data, CLEO [63] 
with a partial- wave analysis of ai — tt~2'k^ decays, and CLEO [65] with a direct partial- 
wave analysis of the r — > UrK^K^n^ decays. The shaded band corresponds to the 
conservative range used in this analysis. 

14.5.3 Total V and A contributions in the non-strange sector 

The total nonstrange vector and axial-vector contributions obtained in this analysis are: 



By = (31.82 ± 0.18 ± 0.12) % , (37) 

Rr,v = 1.787 ±0.011 ±0.007 , (38) 

Ba = (30.19 ± 0.18 ± 0.12) % , (39) 

R^^A = 1.695 ± 0.011 ± 0.007 , (40) 



where the second errors reflect the uncertainties in the V/A separation in the channels 
with KK pairs. Accounting for the correlations between the respective uncertainties, one 
obtains the difference between the vector and axial-vector components, which is physically 
related to the amount of nonperturbative contributions in the nonstrange hadronic r decay 
width (see Section 17): 

Bv-A = (1.63 ± 0.34 ± 0.24) % , (41) 
Rr,v-A = 0.092 ± 0.018 ± 0.014 , (42) 

where again the second error has the same meaning as in Eqs. (37) and (39). The ratio 

p"'^"^ = 0.026 ± 0.007 (43) 
is a measure of the relative importance of nonperturbative contributions. 
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These values yield the normalization of the corresponding spectral functions, which will 
be determined later for the QCD analysis of hadronic r decays and vacuum polarization 
calculations. Compared to the previous analysis, the vector/axial- vector separation is 
significantly improved and will result in better determined nonperturbative contributions. 



14.6 Summary of all measured branching ratios 

A summary list of all r branching ratios measured by ALEPH using the full LEP-1 
statistics is given in Table 17. 



15 Determination of hadronic spectral functions 



15.1 Spectral functions 



The spectral function Vi (ai, ag), where the subscript refers to the spin J of the hadronic 
system, is here defined for a nonstrange vector (axial-vector) hadronic r decay channel 
V~ Vt {A~ Ut). The spectral function is obtained by dividing the normalized invariant 
mass-squared distribution {l/Nv/A){dNv/A/ds) for a given hadronic mass -y/i by the 
appropriate kinematic factor 



vi{s) = 



B{t- V- Vr) 



6 K 



ud\ 



ai{s) = 



B{t- - 

dNy 
Nyds 

B{t- 



1 - 




-1 



(44) 



6|KdP-5Ew B{t- 



X 



ao(s) = 



dNA 
NAds 

B{t- 



1 + 



2s 



mt 



dN, 



A 



6 iKdp 'S'ew B{t e v^v^) Na ds 



mt 



1 - 



mt 



(45) 
(46) 



where ll^^^l = 0.9746 ± 0.0006 [61] denotes the CKM weak mixing matrix element and 
•S'ew = 1-0198 ± 0.0006 accounts for electroweak radiative corrections [50], as discussed 
in Section 16.1. Due to the conserved vector current, there is no J = contribution 
to the vector spectral function, while the only contribution to Oq is assumed to be 
from the pion pole. It is connected via PCAC to the pion decay constant, through 
ao,7r(s) = 47r^/^5(s — m^). The spectral functions are normalized by the ratio of the 
vector/axial- vector branching fraction B{t~ — > V~ /A~ Vt) to the branching fraction of the 
massless leptonic, i.e., electron, channel. For the latter, the value obtained in Section 14.2 
using universality is used, i.e. B{t' — > e" Pe^^r) = (17.810 ± 0.039)%. 
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Table 17: A summary list of ALEPH branching ratios (%). Statistical and systematical 
uncertainties have been combined in quadrature. The labels V, A and S refer to the 
nonstrange vector and axial-vector, and strange components, respectively. CLEO results 
on channels with small branching fractions have been used as indicated. 



mode 

e 

1^ 



17.837 
17.319 



B [%] 
± 0.080 
± 0.077 



TT 




10.828 


± 


0.105 


A 








25.471 


± 


0.129 


V 




71-271^ 




9.239 


± 


0.124 


A 








0.977 


± 


0.090 


V 




7r-47r° 




0.112 


± 


0.051 


A 




7r~7r~7r"'" 




9.041 


± 


0.097 


A 








4.590 


± 


0.086 


V 








0.392 


± 


0.046 


A 








0.013 


± 


0.010 


V 


estimate 


37r-27r+ 




0.072 


± 


0.015 


A 




37r-27r+7r° 




0.014 


± 


0.009 


V 








0.180 


± 


0.045 


V 




(37r)-r/ 




0.039 


d= 


0.007 


A 


CLEO [27] 


«r(^ 7r-7) 




0.040 


± 


0.020 


A 


estimate 


7r'~a;(^ 7r°7, 7r"'"7r' 


-) 


0.253 


± 


0.018 


V 




7r~7r°ci;(-^ 7r°7, tt"^ 


TT-) 


0.048 


± 


0.009 


A 


+ CLEO [26] 


(37r)"cj(— >■ 7r°7, tt" 




0.003 


± 


0.003 


V 


CLEO [27] 






0.163 ± 0.027 


V 








0.145 ± 0.027 


(75 ± 25)% A 








0.153 ± 0.035 


(75 ± 25)% A 








0.163 ± 0.027 


(75 ± 25)% A 




{KKnn)- 




0.050 ± 0.020 


(50 ±50)% A 




K- 
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Using unitarity and analyticity, the spectral functions are connected to the imaginary 
part of the two-point correlation (or hadronic vacuum polarization) functions [67, 68] 

n^^(g) =i/d4^e^^-(0|r(C/^(x)C/^.(0)t)|0) = i-g'^^q' + q'^qn^^iq') +q''q''^%{q') of 
vector (Uj^j = V/^ — qjl^qi) or axial-vector [U^j = Afj = qj^^^^qi) colour-singlet quark 
currents in corresponding quantum states and for time-like momenta-squared q^ > 0. 
Lorentz decomposition is used to separate the correlation function into its J = 1 and 
J = parts. Thus, using the definition (44), one identifies for non-strange quark currents 

Imng^(.) = ^vi{s), 

imn2^(.) = 

ImU^iUs) = ^«o(s), (47) 
which provide the basis for comparing short-distance theory with data. 

15.2 The unfolding procedure 

The measurement of the r spectral functions defined in Eq. (46) requires the determination 
of the invariant mass-squared distributions, obtained from the experimental distributions 
after unfolding from the effects of measurement distortion. The unfolding procedure 
used in this analysis follows a method published in Ref. [69], which has been extensively 
used in previous publications [11, 12]. It is based on the regularized inversion of the 
detector response matrix, obtained from the Monte Carlo simulation, using the Singular 
Value Decomposition (SVD) technique. The regularization function applied minimizes the 
average curvature of the distribution. The optimal choice of the regularization strength 
is found by means of the Monte Carlo simulation where the true distribution is known. 

In order to measure exclusive spectral functions, individual unfolding procedures with 
specific detector response matrices and regularization parameters are applied for each 
r decay channel X considered. An iterative procedure is applied to correct the Monte 
Carlo spectral functions used to subtract the feed across between the modes. Each spectral 
function is determined in 140 mass-squared bins of equal width (0.025 GeV^). 

15.3 Specific systematic studies 

All systematic uncertainties concerning the decay classification are contained in the 
covariance matrix of the branching ratios obtained in Section 10. Therefore only the 
systematic effects affecting the shape of the mass-squared distributions, and not its 
normalization, need to be examined here. 

The general procedure is the same as for the branching ratio analysis. All effects 
affecting the decay classification and the calculation of the hadronic invariant mass are 
considered in turn. Comparisons of data and Monte Carlo distributions are made and the 
corresponding biases are corrected for, as discussed in Section 10, while the uncertainty 
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in the correction is taken as input for the calculation of the systematic uncertainty. In 
the case of the spectral functions, the whole analysis including the unfolding procedure 
is repeated, for each systematic effect. This generates new mass distributions under the 
systematic change which are compared bin-by-bin to the nominal ones, hence providing 
the full 140 X 140 covariance matrix of the spectral function for the studied effect. In 
this process it was found necessary to smooth the mass distributions, before and after 
applying the systematic effect, in order to construct stable covariance matrices. 

Following the systematic studies for the determination of the branching ratios, 
extensive studies are performed to determine the uncertainties at the level of the spectral 
functions. They include the effects from the photon and 7r° energy calibration and 
resolution, the photon detection efficiency (especially in the threshold region above 
350 McV), the shapes of the identification probability distribution, the estimate of the 
number of fake photons, the proximity in the calorimeter of other photon showers and 
of energy deposition by charged particles, and the separation between radiative and 7r° 
decay photons for residual single photons. The photon and 7r° reconstruction constitutes 
the main part of the uncertainty for the determination of spectral functions. 

Similarly, the effects of momentum calibration and resolution uncertainties in the 
reconstruction of charged tracks are checked, accompanied by tests of the reconstruction 
efficiency of highly coUimated multi-prong events, and the simulation of secondary nuclear 
interactions. 

In addition, systematic errors introduced by the unfolding procedure are tested by 
comparing known, true distributions to their corresponding unfolded ones. 

Finally, systematic errors due to the limited Monte Carlo statistics and to uncertainties 
in the branching ratios are added. 

15.4 Spectral functions for nonstrange exclusive modes 
15.4.1 The results 

Before unfolding the mass distributions, the r and non-r backgrounds are subtracted. 
In the case of r feedthrough the Monte Carlo distributions normalized to the measured 
branching fractions are used. Since the spectral functions are measured for nonstrange 
exclusive final states, the contributions from strange modes classified in the same topology 
are subtracted using their Monte Carlo spectral functions normalized by the measured 
branching fractions [9]. The measurement of the strange spectral functions has been 
pubhshed elsewhere [9]. 

The vector spectral functions of the dominant two- and four-pion modes are given in 
Figs. 54, 55 and 56. The error bars shown correspond to the diagonal elements of the 
covariance matrix. They include both statistical and systematic uncertainties. It should 
be pointed out that the unfolding generates strong correlations between neighbouring 
bins. This information is contained in the covariance matrix that is determined for each 
spectral function ^. 

^The corrected invariant mass-squared distributions from this analysis and their covariance matrices 
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Figure 54: The spectral function for the tttt^ hadronic mode. 

Figures 57 and 58 show the unfolded 2tt~tt~^ and 7r~27r° z/^ spectral functions which 
are the dominant contributions to the axial-vector spectral function. Both spectra 
are dominated by the ai resonance and indeed the two distributions have consistent 
shape, while their normalization was already observed to agree, as expected from isospin 
invariance with a small correction from the slightly different charged and neutral pion 
masses. The comparison is seen in Fig. 59. 

The 5-pion channels, n'iiT^i'T-, 27r~7r"'"27r°i/T- and 37r~27r+z/T-, have small branching 
ratios and relatively large feedthrough background (except for the 5-charged-pion mode). 
Their combined mass-squared distribution, shown in Fig. 60, does not agree with 
the TAUOLA decay simulation in the KORALZ07 generator [18]. The experimental 
spectrum is harder than the phase-space model used in TAUOLA, in the absence of prior 
experimental information on the dynamics in these decay modes. 

15.4.2 Comparison with other experiments 

Spectral functions for some specific final states have been determined by CLEO and 
OPAL. Since the data are not presented with the same binning and furthermore plotted 
sometimes against the mass rather than the mass-squared, the comparison must involve 
some treatment of the original data. 

are publicly accessible at the website http://aleph.web.lal.in2p3.fr/tau/specfun for the inclusive as well 
as the main exclusive channels. 
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Figure 55: The spectral function for the Stttt'^ hadronic mode. While the original spectral 
function is determined in 0.025 GeV^ bins, it is rebinned here to 0.1 GeV^ for easier 
reading. 
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Figure 56: The spectral function for the ttStt^ hadronic mode. While the original spectral 
function is determined in 0.025 GeV^ bins, it is rebinned here to 0.1 GeV^ for easier 
reading. 
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Figure 57: The spectral function for the Stt hadronic mode. 
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Figure 58: The spectral function for the 7r27r° hadronic mode. 
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Figure 59: Comparison of the corrected mass spectra for the Svr and 7r27r° hadronic modes. 

The results for the 7r7r° spectral function from ALEPH, CLEO [58], and OPAL [70] 
are compared in Fig. 61 in the p resonance region. For this comparison, each data 
set is normalized to the weighted mean of the three measurements. The different 
data are in agreement within their quoted uncertainties, when taking into account the 
correlation between different mass bins, particularly the strong statistical correlation 
between neighbouring points due to the unfolding procedure. It should be pointed out that 
each spectral function is normalized to the world-average branching ratio for r — > v^tttt^ 
that is dominated by the ALEPH result. Hence the present exercise should be understood 
as a test of the shape of the spectral function. However, if the branching ratio determined 
by each experiment (see Section 16.2.3) is used instead, the spectral functions are still in 
agreement, albeit with larger uncertainties in the case of CLEO and OPAL. The two most 
precise results from ALEPH and CLEO do agree well. The statistics is comparable in 
the two cases, however due to a fiat acceptance in ALEPH and a strongly increasing one 
in CLEO, ALEPH data are more precise below the p peak, while CLEO is more precise 
above. 

15.5 Results on inclusive V ^ A and V :^ A spectral functions 
15.5.1 The vector spectral function 

The inclusive r vector spectral function with its most important contributions is shown 
in Fig. 62. The dashed line depicts the naive parton model prediction while the massless 
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Figure 60: Comparison of the background-subtracted mass-squared spectrum in the 5- 
pion modes {h^Ah^Ut-, 27r^7r+27r°z/^ and 37r^27r^z/^) with the phase-space model used in 
the TAUOLA decay generator [18]. The dashed histogram is the generated distribution, 
while the solid histogram is the reconstructed spectrum with the detector response which 
can be directly compared to the data (points). A clear disagreement is observed. 



102 




_ 0.3 
(-> 

0.2 

1 0.1 
CO 

0) 

cp 

\1 

u_'' -0.1 
I 

E -0.2 

LL 

^ -0.3 



0.5 



ALEPH 

CLEO 

OPAL 




0.55 



0.6 



s (GeV^) 



0.65 



0.7 



Figure 61: Relative comparison of the tt+tt^ spectral functions extracted from r data 
from different experiments, expressed as a ratio to the average r spectral function. The 
lower figure emphasizes the p region. For CLEO only statistical errors are shown. 
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Figure 62: The total vector spectral function. The shaded areas indicate the contributions 
from the exclusive r vector channels, where the shapes of the contributions labelled 'MC 
are taken from the Monte Carlo simulation. The lines show the predictions from the naive 
parton model and from massless perturbative QCD using a;s(M|) = 0.120. 
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perturbative QCD prediction [71] using as^M^) = 0.120 (solid line) lies roughly 25% 
lower than the data at m^. Although the statistical power of the data is weak near the 
kinematic limit, the trend of the spectral function clearly indicates that the asymptotic 
region is not reached. 

The two- and four-pion final states are measured exclusively, while the six-pion state 
is only partly measured. The total six-pion branching ratio has been determined in [11] 
using isospin symmetry. However, one has to account for the fact that the six-pion channel 
is contaminated by isospin- violating t~ ^ r]2-K~T:^ r)7c~2'K^i'r decays, as reported by 
the CLEO Collaboration [27]. 

The small fraction of the uj Vj. decay channel that is not reconstructed in the four- 
pion final state is added using the simulation. Similarly, one corrects for rjTi^n^ Vt decay 
modes where 77 decays into pions. For the 77 — > 27 mode, the r decay is classified in 
the /iSTT^i^T- final state, since the two-photon mass is inconsistent with the 7r° mass and 
consequently each photon is reconstructed as a 7r°. 

The K~ Vt mass distribution is taken entirely from the simulation. The vector 
contributions in the KKti and KKnn modes are taken from Section 14.5. The 
corresponding spectral functions are obtained from the Monte Carlo simulation. 

The invariant mass spectra of the small contributions labelled 'MC in Figs. 62 and 
63 are taken from the Monte Carlo simulation accompanied by a channel-dependent 
systematic error of up to 50% of the bin entry. 



15.5.2 The cLxial-vector spectral function 

In complete analogy to the vector spectral function, the inclusive axial-vector spectral 
function is obtained by summing up the exclusive axial- vector spectral functions with the 
addition of small unmeasured modes taken from the Monte Carlo simulation. 

The small fraction of the a;7r~7r° 1/^ decay channel that is not accounted for in the 
27r'~7r''"27r° u^- final state is added from the simulation. Also considered are the axial- 
vector ?7(37r)^z/^ final states [27]. CLEO observed that the dominant part of this 
mode issues from the ^fi(1285)7r~ intermediate state, with B{t~ — s>fi7r~z/^) = 
(0.068 ± 0.030)%, measured in the fi — >• t^tt+tt" and fi — >• t^tt^tt^ decay modes [27]. 
Since the fi meson is isoscalar, the branching ratios relate as S(t~ r]2T:~T:^ v^-) — 
2 X Bir" — > r77r~27r° i^T-). The distributions are taken from the ordinary six-pion phase 
space simulation accompanied by large systematic errors. As discussed for the vector 
spectral function, the KKti final states contribute dominantly to the inclusive axial- 
vector spectral function, with full anticorrelation to the inclusive vector spectral function. 
Invariant mass distributions for these channels are taken from the simulation. 

The total inclusive axial-vector spectral function is plotted in Fig. 63 together with 
the naive parton model and the massless, perturbative QCD prediction. One observes 
that the asymptotic region is apparently not reached at the r mass scale. 
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Figure 63: The total axial-vector spectral function. The shaded areas indicate the 
contributions from the exclusive r vector channels, where the shapes of the contributions 
labelled 'MC are taken from the Monte Carlo simulation. The lines show the predictions 
from the naive parton model and from massless perturbative QCD using ^^(Ml) = 0.120. 
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Figure 64: The inclusive vector plus axial- vector [vi+ai) spectral function and predictions 
from the parton model and from massless perturbative QCD using as(M|) = 0.120. 

15.5.3 The V ± A spectral functions 

For the total {vi + ai) hadronic spectral function one does not have to distinguish the 
current properties of the non-strange hadronic r decay channels. Hence the mixture of all 
contributing non-strange final states is measured inclusively using the following procedure. 

The two- and three-pion final states dominate and their exclusive measurements are 
added with proper accounting for the correlations. The remaining contributing topologies 
are treated inclusively, i.e., without separation of the vector and axial- vector decay modes. 
This reduces the statistical uncertainty. The effect of the feedthrough between r final 
states on the invariant mass spectrum is described by the Monte Carlo simulation and 
thus corrected in the data unfolding. In this procedure the simulated mass distributions 
are iteratively corrected using the exclusive vector /axial- vector unfolded mass spectra. 
Another advantage of the inclusive {vi + ai) measurement is that one does not have to 
separate the vector /axial- vector currents of the KKvr and KKvrvr modes. The {vi + ai) 
spectral function is depicted in Fig. 64. The improvement in precision when comparing 
to a sum of the two parts (Fig. 62 and Fig. 63) is obvious at higher mass-squared. 

One clearly sees the oscillating behaviour of the spectral function but, unlike the 
vector/axial- vector spectral functions, it does approximately reach the asymptotic limit 
predicted by perturbative QCD at s — >• m^. Also, the V + A spectral function, 
including the 5-function vr contribution, exhibits the features expected from global quark- 
hadron duality: despite the huge oscillatory effects from vr, and p, ai and p' hadron 
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Figure 65: The vector minus axial- vector (f i — ai) spectral function. In the parton model 
as well as in perturbative QCD vector and axial- vector contributions are degenerate. 

resonances, the spectral function qualitatively averages out to the quark contribution 
from perturbative [72]. This observation is the physical basis for the quantitative QCD 
analysis performed in Section 17. 

In the case of the {vi — Oi) spectral function, uncertainties on the V/A separation 
are reinforced due to their anticorrelation. In addition, anticorrelations in the branching 
ratios between r final states with adjacent numbers of pions increase the errors. The 
{vi —ai) spectral function is shown in Fig. 65. The oscillating behaviour of the respective 
Vi and fli spectral functions is emphasized and the asymptotic behaviour is clearly not 
reached at m^. However again here, the oscillation generated by the hadron resonances 
qualitatively averages out to zero, as predicted by perturbative QCD. 



16 The vector spectral functions and tests of CVC 

16.1 General remarks on spectral functions in r decays and e+e~ 
annihilation 

Using the optical theorem and in the limit of isospin invariance (implicit in the hypothesis 
of CVC — Conserved Vector Current) the spectral function of a r decay mode X~ i/-,-, where 
the hadronic system X is in a vector state, is related to the e'^e~ annihilation cross section 
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of the corresponding isovector final state X": 

ai+l-^xo = ^^i,x- , (48) 

with the electromagnetic fine structure constant a. 

Since the breaking of isospin symmetry is expected at some level, in particular from 
electromagnetic effects, it is useful to carefully write down all the factors involved in the 

comparison of e~^e~ and r spectral functions in order to make explicit the possible sources 
of CVC breakdown. For the dominant tttt spectral functions, one has on the e+e" side 

(T(e+e- — ^TT+TT") = ^^vo{s) , (49) 
-o(.) = ^\F'As)\', 

where /3q(.s) is the threshold kinematic factor and F^{s) the pion form factor. On the r 
side, the physics is contained in the hadronic mass distribution through 



1 dr 

f dj' 



— — IT — > 71 n^u-r 



-C{s)v-{s) , (50) 



s V f 2s 



ml I \ ml 



SU(2) symmetry imphes v^{s) — Vq{s). The threshold functions /3o,- are defined by 

/3o,- = /5(s, m^- , m^o,- ) , (51) 

where 



/3(s,mi,m2) 



' ^ _ {mi + m2Y \ L _ (mi - 7712)^ ' 



1/2 



(52) 



In this analysis of r decays the rate is measured inclusively with respect to radiative 
photons, i.e., for r z/^7r7r°(7). The measured spectral function is thus v*_{s) — 
V-{s) G{s), where G{s) is a radiative correction. 

Several levels of SU(2) breaking can be identified. 



• Electroweak radiative corrections to t decays are contained in the 5'ew factor [50, 73] 
which is dominated by short- distance effects. As such it is expected to be 
weakly dependent on the specific hadronic final state, as verified in the r — > 
{t:,K)vt channels [53]. Recently, detailed calculations have been performed for 
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the TTTT^ channel [24], which also confirm the relative smallness of the long-distance 
contributions. The total correction is 

chad Chad 

^EW - -ji^ [06) 

'^EM 

where S^^r is the leading-log short- distance electroweak factor (which vanishes for 
leptons) and S'^'^'^'^^ are the nonleading electromagnetic corrections. The latter 
corrections are calculated in Ref. [73] at the quark level and in Ref. [24] at the 
hadron level for the tttt^ decay mode, and in Refs. [50, 73] for leptons. The total 
correction amounts [57] to 5"^^" = 1.0198 ± 0.0006 for the inclusive hadron decay 
rate and = (1.0232 ± 0.0006) G][^{s) for the nn^ decay mode, where G'l;';^(s) 
is an s-dependent radiative correction [24]. 

• The pion mass splitting breaks isospin symmetry in the spectral functions [54, 74] 
since P-{s) ^ Po{s)- 

• Isospin symmetry is also broken in the pion form factor [54, 24] from the tt mass 
splitting. 

• A similar effect is expected from the p mass splitting. The theoretical 
expectation [75] gives a limit (< 0.7 MeV), but this is only a rough estimate. 
Hence the question must be investigated experimentally, the best approach being 
the explicit comparison of r and e^e~ 27r spectral functions, after correction for the 
other isospin-breaking effects. No correction for p mass splitting is applied initially. 

• Explicit electromagnetic decays such as 777, 777, l'^l~ and 77777 introduce small 
differences between the widths of the charged and neutral p's. 

• Isospin violation in the strong amplitude through the mass difference between u and 
d quarks is expected to be negligible. 

• An obvious, but large correction must be applied to the r spectral function to 
introduce the effect of p — a; mixing, only present in the neutral channel. This 
correction is computed using the parameters determined in the e'^e~ experiments 
in their fits of the 7r"'"7r~ lineshape taking into account the p — cu interference [60]. 



16.2 Comparison with e'^e~ results 
16.2.1 The 27r spectral function 

Figure 66 shows the comparison for the 27r spectral functions obtained by averaging the 
results of the present analysis and the published data from CLEO [58] and OPAL [70]. 
The e+e" data are taken from TOE [76], OLYA [77], CMD [78], CMD-2 [60], DM1 [79], 
DM2 [80], and NA7 [81]. The most precise results come from CMD-2 and have undergone 
recently a major revision due to previously incorrect radiative corrections [59]. Visually, 
the agreement seems satisfactory, however the large dynamical range involved does not 
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permit an accurate test. To do so, the e+e~ data are plotted as a point-by-point ratio to 
the T spectral function in Fig. 67, also showing the recent accurate data from KLOE [82] 
using the radiative return technique. The data are found to be consistent below and 
around the p peak, while a discrepancy is observed for masses larger than 0.85 GeV. 



16.2.2 The 47r spectral functions 

The spectral function measurements of the r vector current final states 7r~37r° and 
27r"7r"'"7r° are compared to the cross sections of the corresponding e+e" annihilation into 
the isovector states 2t^~2t^^ and 7r~7r+27r°. Using Eq. (46) and isospin invariance the 
following relations hold: 

(Je+e--+7r+7r-7r+7r- = 2- — V^^^^- ^^^Q „^ , (54j 

i=i - ^^"^ r 1 

s 

The comparison of the Att cross sections is given in Fig. 68 for the 2tt~^2tt~ channel 
and in Fig. 69 for 7r"'"7r~27r'', where the latter suffers from large differences between 
the results from the various e~^e~ experiments. The r data, combining two measured 
spectral functions according to Eq. (55) and corrected for isospin breaking originating 
from the charged- vs- neutral pion mass difference [74] , lie somewhat in between with large 
uncertainties above 2 GeV^ because of the lack of statistics and a large feedthrough 
background in the r VrTT~^Ti^ mode. In spite of these difficulties the tt^Stt" spectral 
function is in agreement with e^e~ data as can be seen in Fig. 68. It is clear that intrinsic 
discrepancies exist among the e^e~ experiments and that a quantitative test of CVC in 
the 7r'^7r~27r° channel is premature. 



16.2.3 Branching Ratios in r Decays and CVC 

It is possible to compare the e^e~ and r spectral functions in a more quantitative way 
by calculating weighted integrals over the mass range of interest up to the r mass. One 
convenient choice is provided by the r branching fractions which involve as a weight the 
kinematic factor (1 — ^)^(1 -|- ^) coming from the V — A charged current in r decay. It 
is then possible to directly compare the measured r branching ratios to their prediction 
through isospin invariance (CVC) using as input the e"'"e~ isovector spectral functions. 

Using the "universality-improved" branching fraction given in Eq. (32), the results for 
the main channels are given in Table 18. The errors quoted for the CVC values are split 
into uncertainties from [i) the experimental input (the e^e~ annihilation cross sections) 
and the numerical integration procedure, {ii) additional radiative corrections applied to 
some of the e+e" data (see Ref. [57] for details), and {iii) the isospin-breaking corrections 
when relating r and e+e" spectral functions. 

As expected from the preceding discussion, a large discrepancy is observed for the r — > 
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Figure 66: Comparison of the tt+tt^ spectral functions from e^e~ and isospin-breaking 
corrected r data from ALEPH, CLEO, and OPAL expressed as e^e~ cross sections. The 
band indicates the combined e^e~ and r result within la errors. It is given for illustration 
purpose only. The complete references for the e^e~ data are given in Refs. [57, 61] and 
[81]. 
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Figure 67: Relative comparison of the n~^7r~ spectral functions from e~^e~ and isospin- 
breaking corrected r data from ALEPH, CLEO, and OPAL expressed as a ratio to the r 
spectral function. The band shows the uncertainty on the latter function. The complete 
references for the e~^e~ data are given in Refs. [57, 61] and [82]. 
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Figure 68: Comparison of the 27r+27r^ spectral functions from e~^e~ and isospin-breaking 
corrected r data, expressed as e"'"e~ cross sections. The complete references for the e~^e~ 
data are given in Refs. [57, 61]. 
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Figure 69: Comparison of the 7r^7r^27r spectral functions from e^e~ and isospin-breaking 
corrected r data, expressed as e^e^ cross sections. The complete references for the e~^e~ 
data are given in Refs. [57, 61]. 
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Table 18: Branching fractions of r vector decays into 2 and 4 pions in the final state. 
Second column: present ALEPH results. Third column: inferred from e^e~ spectral 
functions using the isospin relations (48,50,51,54, 55) and correcting for isospin breaking 
following Refs. [57, 61]. Experimental errors, including uncertainties on the integration 
procedure, and theoretical (missing radiative corrections for e^e~ , and isospin-breaking 
corrections and Vud for r) are shown separately. Right column: differences between the 
direct measurements in r decays and the CVC evaluations, where the separate errors have 
been added in quadrature. 





Mode 


T ALEPH 


Branching fractions (in %) 
e+e- via CVC 


A(r- e+e~) 




" — ^ V-,-1T~TT^ 


25.47 ±0.13 


24.52 ± 0.26exp ± O.llrad ± 0.12su(2) 


+0.95 ± 0.33 








^ V ' 

0.31 




T~ 


' — ^ VtTT^^'K^ 


0.98 ± 0.09 


1.09 ± O.OGexp ± 0.02rad ± 0.05su(2) 


-0.11 ±0.12 








0.08 




T' 




4.59 ± 0.09 


3.63 ± 0.19,,xp ± 0.0 ± 0.09su(2) 


+0.96 ± 0.23 








0.21 





Vr 7r~7r° branching ratio, with a difference of (0.95 ± 0.13^ ± 0.26ee ± O.llrad i 0.12su(2))%, 
where the uncertainties are from the r branching ratio, e^e~ cross sections, e^e~ missing 
radiative corrections and isospin-breaking corrections (including the uncertainty on Vud)-, 
respectively. Adding all errors in quadrature, the effect represents a 2.9 a discrepancy. 
Since the disagreement between e^e" and r spectral functions is more pronounced at 
energies above 0.85 GeV, a smaller discrepancy is expected in the calculation of the 
hadronic contribution to the muon anomalous magnetic moment because, in this case, 
the spectral function is weighted by a steeply faUing kernel K{s). 

Besides the result from the present analysis, a similar comparison can be performed 
using results from other experiments. The values for the r branching ratios involve 
measurements [37, 36] given without charged hadron identification, i.e. for the hT^^Uri 
h^TT^UT- and 3/i7r*^z/^ final states. The corresponding channels with charged kaons have 
been measured [9, 83] and their contributions are subtracted out in order to obtain the 
pure pionic modes. The comparison is displayed in Fig. 70 where it is clear that ALEPH, 
CLEO, and OPAL all separately disagree with the e+e~-based CVC result. 

The situation in the An channels is different. Agreement is observed for the 7r~37r° 
mode within a relative accuracy of 12%, however the comparison is not satisfactory for 
the 2n~7i~^7i^ mode. In the latter case, the relative difference is very large, (23 ± 6)%, 
compared to a reasonable level of isospin symmetry breaking. As such, it rather points 
to experimental problems that have to be investigated, which are emphasized by the 
scatter observed among the different e'^e~ results. In this respect, it is to be noted 
that new preliminary results from the SND and CMD-2 experiments have been recently 
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e"'e CVC I — • — I ^ average 
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23 24 25 26 27 
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Figure 70: The measured branching ratios for r —>■ z/^7r~7r° compared to the prediction 
from the e~^e~ tt^tt^ spectral function applying the isospin-breaking correction factors 
discussed in Section 16.1. The measured branching ratios are from ALEPH, CLEO [37] 
and OPAL [36]. The OPAL result was obtained from their /ivr" branching ratio, reduced 
by the small Ktt^ contribution measured by ALEPH [9] and CLEO [83]. 

presented [84]; while the new data contradict earlier results from the same experiments, 
they are in good agreement with the ALEPH spectral functions. 



16.3 Fits to the tttt^ spectral function 
16.3.1 Procedure 

The TTTT^ spectral function is dominated by the wide p resonance, parametrized in the 
forthcoming fits following Gounaris-Sakurai [85] (GS). The GS parametrization takes into 
accounts analyticity and unitarity properties. The fits make use of the covariance matrix 
taking into account the correlations between different mass bins. 

If one assumes vector dominance, the pion form factor is given by interfering 
amplitudes from the known isovector meson resonances p(770), /9(1450) and p(1700) with 
relative strengths 1, /?, and 7. Although one could expect from the quark model that /3 
and 7 are real and respectively negative and positive, the phase of (3, (pp is left free in 
the fits, while the much smaller 7 is assumed to be real for lack of precise experimental 
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information at large masses. Taking into account p — uj interference, one writes 

l+5BW^(7S3)(s) 



r.i'=l,0. 



s) 



BW 



p(770) 



with the Breit-Wigner propagators 



where 



1 + /3 + 7 



ml{l + d-Tp/mp) 



(56) 



m? - s + f{s) - i^/s^p{s) ' 



(57) 



e{s) {h{s)-h{ml)) + {ml-s)k'{ml) 



The P-wave energy-dependent width is given by 



dh 
ds 



(58) 



Tp{s) = Tp{ml) 



(59) 



where k{s) = | \/s P (s) and A;(mp are pion momenta in the p rest frame. The function 
h{s) is defined as 



^ 2 *Wi„yi±2M£) 

TT ^/S 2m7r 



(60) 



with dh/ds\ 



h{ml) (8A;2(m2))-i - (2m2) 



2^-1 



{2-Km ) . Interference with the 



isospin-violating electromagnetic uj tt+tt^ decay occurs only in e~^e~ annihilation. 
Consequently, 6 is fixed to zero when fitting r data. The normalization BW^(^^)(0) = 1 
fixes the parameter d — f{0)/{Tpmp), which is found to be [85] 



3 ml ^^mp + 2k{ml) ^ 
7rk^{ml) 2m7r 



+ 



m^mp 



2tt k{mV) 7rk^{mV) 



(61) 



16.3.2 Fit to the ALEPH data 

The result of the fit to the ALEPH data using the GS paramctrization is given in Table 19 
and illustrated in Fig 71. The p mass uncertainty is dominated by systematic effects, the 
largest being the knowledge of the 7r° energy scale (calibration). 

Concluding from Table 19, the fit establishes a need for the p(1450) contribution to 
the weak pion form factor. Weak evidence is found for a p(1700) contribution which is 
expected since the resonance hes close to the spectrum end-point. Most of the fitted 
parameters exhibit large correlations as seen in Table 20. 
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Table 19: Fit results of the ALEPH pion form factor in t~ 7r~7r° z/^ decays using the 
Gounaris-Sakurai (GS) parametrization. The parameters rnp±(^noci) and rp(i7oo) are kept 
fixed to values obtained from fits of e~^e~ data extending in mass-squared up to 3.6 GeV^. 



Paranioter 


ALEPH GS 


"^p±(770) (MeV) 
r^±(770) (MeV) 


775.5 ± 0.7 
149.0 ± 1.2 


P 

(f)B (degrees) 
^^ip±(i450) (MeV) 
rp(i450)(MeV) 


0.120 ± 0.008 
153 ± 7 
1328 ± 15 

468 ± 41 


7 

"^P±(i700) (MeV) 
Tpdroo) (MeV) 


0.023 ± 0.008 
[1713] 
[235] 


XVDF 


119/110 



Table 20: The correlation coefficients between the fitted parameters in the fit of the 
ALEPH 7r7r° spectral function. 



ALEPH 


"^p(770) 


rp(770) 




4>I3 


^p(1450) 


rp(1450) 


7 


"^p(770) 


1 


0.38 


0.29 


-0.28 


0.42 


0.46 


-0.36 


rp(770) 




1 


0.52 


-0.02 


0.16 


0.71 


-0.38 


P 






1 


0.31 


-0.08 


0.63 


0.03 


0/3 








1 


-0.87 


-0.40 


0.77 


^p{1450) 










1 


0.58 


-0.94 


rp(1450) 












1 


-0.66 


7 














1 



16.4 Combined fit to r and e+e data 

It is interesting to perform a combined fit using r (ALEPH and CLEO) and e'^e~ data in 
order to better constrain the lesser known parameters in the phenomenological form factor. 
For this study, the r spectral function is duly corrected for the isospin-breaking effects 
identified in Section 16.1. Also photon vacuum polarization contributions are removed in 
the e~^e^ spectral function since they are absent in the r data. In this way, the mass and 
width of the dominant p(770) resonance in the two isospin states can be determined. For 
the sub-leading amplitudes from the higher vector mesons, isospin symmetry is assumed 
and therefore common masses and widths are used in the fit. 

The result of the combined fit is given in Table 21. The differences between the masses 
and widths of the charged and neutral p(770)'s can be extracted to yield 

rup- - rupo = (2.4 ± 0.8) MeV (62) 



119 




Figure 71: Fit of the ALEPH tt tt*^ form factor squared using the Gounaris-Sakurai 
parametrization. 

r^. - Tpo = (0.2 ± 1.0) MeV (63) 

The mass sphtting is somewhat larger than the theoretical prediction (< 0.7 MeV) [75], 
but only at the 2a level. The expected width splitting, from known isospin breaking, but 
not taking into account any p mass splitting, is (0.7 ± 0.3) MeV [24, 57]. However, if the 
mass difference is taken as an experimental fact, then a larger width difference would be 
expected. From the chiral model of the p resonance [86, 24], one expects 

where Ar^;^ is the width difference from electromagnetic decays (as discussed above), 
leading to a total width difference of (2.1 ± 0.5) MeV, marginally consistent with the 
observed value. 

Mass splitting for the p was in fact considered in our previous analysis [3]. Using 
pre-CMD-2 e~^e~ data, a combined fit was attempted which produced a mass splitting 
consitent with within an uncertainty of 1.1 MeV. However, the form factor from e~^e~ 
data still contained the vacuum polarization contribution (producing a 1.1 MeV shift) and 
a normalization problem was discovered in the treatment of the r data in the combined 
fit. With the advent of precise CMD-2 data [59], it became apparent that differences 
were showing up between r and e~^e~ form factors. A large part of the discrepancy was 
removed when CMD-2 re-analysed their data [60]. Since the r results from ALEPH, 
CLEO and OPAL have been shown to be consistent, and since the recent results from 
the radiative return analysis of KLOE [82] are in fair agreement with the corrected 
CMD-2 results, the question of the p mass splitting is now more relevant. However a 
combined fit as shown in Table 22 using the r data on one side and the e~^e~ data on 
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Table 21: Combined fit to the pion form factor-squared to ALEPH, CLEO r and 
all e^e~ data, where vacuum polarization has been excluded in the latter data. The 
parametrization of the p(770), p(1450), p(1700) line shapes follows the Gounaris-Sakurai 
formula. Separate masses and widths are fitted for the p(770), while common values are 
kept for the higher vector mesons. All mass and width values are in units of MeV and 
the phases are in degrees. 





T and e^e 


"^p-{770) 


1 i o.o in u.u 




773.1 ± 0.5 


rp-(770) 


148.2 ± 0.8 


rp()(770) 


148.0 ± 0.9 




(2.03 ± 0.10) 10-^ 




(13.0 ±2.3) 


P 


0.166 ± 0.005 




177.8 ± 5.2 


"^'p(1150) 


1409 ± 12 


^^(1150) 


501 ± 37 


7 


0.071 ± 0.006 


4>i 


[0] 


"T-p(1700) 


1740 ± 20 


rp(1700) 


[235] 


XVDF 


383/326 



the other, and requiring for consistency the constraint from Eq. (64) has a probability 
of only 0.6%. In fact, correcting for different masses extracted from the fit and using 
the corresponding constrained widths, improves the agreement between the r and e+e~ 
line shapes, but at the expense of a significant discrepancy in normalization. It should be 
noted that an additional correction for this apparent p mass splitting increases the present 
discrepancy for the muon anomalous magnetic moment between the estimates based on 
r and e^e~ spectral functions [87]. Looking back at the detailed comparison between 
spectral functions in Fig. 67, it is seen that the new KLOE data shows a trend below 
and above the p resonance which is not in good agreement with CMD-2 (and even worse 
with T data). This behaviour prevents one at the moment from attributing the e'^e'-T 
difference to a simple isospin-breaking correction resulting from different p masses and 
widths. The origin of the discrepancy thus remains ununderstood. 



17 QCD analysis of r spectral functions 

The measurement of as{m'^) presented in this section is an update with the full r sample 
of the analysis pubhshed in Ref. [12]. The method is based on a simultaneous fit of 
QCD parametrizations including perturbative and nonperturbative components to the 
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Table 22: Combined fit to the pion form factor squared to ALEPH, CLEO r and 
e+e" data, where vacuum polarization has been excluded for the latter data. The 
paramctrization of the p(770), p(1450), p(1700) line shapes follows the Gounaris-Sakurai 
formula. For the p(770), only nip-, nipO, and Fp- are fitted, while FpO is computed from 
Eq. (64). All mass and width values are in units of MeV and the phases are in degrees. 
The parameters related to p(1450) and p(1700) amplitudes are fitted, assuming they are 
identical in both data sets. 





T and e^f 


'Tip- (770) 


775.4 ± 0.6 


"ipO(770) 


773.1 ± 0.5 


rp-(770) 


148.8 ± 0.8 


FpO (770) 


(146.7) 


QJpoj 


(2.02 ± 0.10) 10-^ 




(15.3 ±2.0) 


P 


0.167 ± 0.006 


4>f3 


177.5 ± 6.0 


"ip(1450) 


1410 ± 16 


rp(1450) 


505 ± 53 


7 


0.071 ± 0.007 




[0] 


W-p(1700) 


1748 ± 21 


Fp(1700) 


[235] 




390/327 
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measured ratio Rt 

= r(r-^hadrons-..) 

and to the spectral moments defined below (Section 17.2). It was proposed by Le Diberder 
and Pich [88] and has been employed in previous analyses by the ALEPH [10, 12], 
CLEO [89] and OPAL [70] Collaborations. 



17.1 Theoretical prediction 

According to Eq. (47) the absorptive parts of the vector and axial-vector two-point 
correlation functions n^'^^y^^(s), with the spin J of the hadronic system, are proportional 
to the T hadronic spectral functions with corresponding quantum numbers. The non- 
strange ratio Rt- can be written as an integral of these spectral functions over the invariant 
mass-squared s of the final state hadrons [67] : 



Rr = 127r5Ew 1 + Imn(i)(s + le) + Imn(°)(s + 

J mi;. \ J [ \ J 



ie) 



(66) 



where H^'^^ can be decomposed as li^'^'' — {p-uly + '^^uI^a)- The correlation function 

n'-'^^ is analytic in the complex s plane everywhere except on the positive real axis where 
singularities exist. Hence by Cauchy's theorem, the imaginary part of 11*^ "'^ is proportional 
to the discontinuity across the positive real axis. 

The energy scale Sq = ""^r large enough that contributions from nonpcrturbative 
effects are small. It is therefore assumed that one can use the Operator Product Expansion 
(OPE) to organize perturbative and nonpcrturbative contributions to i?r(so). The factor 
(1 — s/rri^Y suppresses the contribution from the region near the positive real axis where 
n("^)(s) has a branch cut and the OPE validity is restricted [90]. 

The theoretical prediction of the vector and axial-vector ratio Rt,v/a can thus be 
written as: 

Rry,A = ||K.r^Ew fl + + + + E €y/^ > (67) 

^ V D=4,6,... / 

with the residual non-logarithmic electroweak correction (5^w = 0.0010 [73], neglected in 
the following, and the dimension D — 2 contribution S^^ vva from quark masses which is 
lower than 0.1% for u, d quarks. The term S^^^ is the massless perturbative contribution, 
while the S^^^ are the OPE terms in powers of : 

xiD) _ \- r AOud{^^))v/A , . 

^udyiA - 1^ <^udy/A[s,l^) , ' (68) 

where the parameter /i separates the long-distance nonpcrturbative effects, absorbed into 
the vacuum expectation elements (C„(i(/x)), from the short-distance effects which are 
included in the Wilson coefficients Cudy/A{s, /jl) [91]. 
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The perturbative prediction adopted in this analysis follows in detail Ref. [92]. The 
perturbative contribution is given in the chiral limit. Effects from quark masses have 
been calculated in Ref. [93] and are found to be well below 1% for the light quarks. Thus 
the contributions from vector and axial-vector currents coincide to any given order of 
perturbation theory and the results are flavour independent. 

The perturbative contribution in Eq. (67) is then given by [92] 



1 + 5(0) ^ , 



(69) 



n=0 



with Kq = Ki = 1, K2 = 1.63982 and K3 = 6.37101 for three active flavours in the MS 
scheme [71]. The coefficients Kn arc known up to three-loop order and for n > 2 they 
depend on the renormalization scheme employed. Although progress was recently made 
in the evaluation of some contributions at the four-loop order [94] , the full calculation of 
K4 is not yet available. The functions A^") (ag) in Eq. (69) are the contour integrals 



^ 2m J s 

\s\=ml 




TT 



(70) 



where the contour runs counter clockwise around the circle from -|- ie to — ie. The 
strong coupling constant in the vicinity of can be expanded in powers of q;s(so), with 
coefficients that are polynomials in ln(s/m^) [67]. The perturbative prediction becomes 
then a function of the coefficients and elementary integrals. Up to fourth order the 
fixed-order perturbation theory (FOPT) expansion reads 



1 + 5^^^ = 1 + ^^ii!^ + 5.2023 f ^iMlY + 26.366 ("^i^T 



TT 



+ (i^4 + 78.00) 



TT 



TT 



TT 



(71) 



with the unknown Kt^ coefficient. 

Another approach to the solution of the contour integral (70) is to perform a direct 
numerical evaluation using the solution of the renormalization group equation (RGE) 
to four loops [98] as input for the running as{—s) [99, 92]. It implicitly provides a 
resummation of all known higher order logarithmic integrals and improves the convergence 
of the perturbative series. While, for instance, the third order term in the expansion (71) 
contributes with 17% to the total (truncated) perturbative prediction, the corresponding 
term of the numerical solution amounts only to 6.6% (assuming ^^(m^) = 0.35). 
This numerical solution of Eq. (69) will be referred as contour-improved fixed-order 
perturbation theory (FOPTci) in the following. 



Despite a number of arguments expressed in Ref. [92] , the intrinsic ambiguity between 
FOPT and FOPTci is unresolvable at present. This is due to the truncation of the 
perturbative approximation of 5^^^ at finite order in ctg. A conservative measure of this 
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ambiguity is obtained from the deviation in found when cutting all additional orders 
in as (which is FOPT) and keeping them (FOPTci), respectively. Both methods are 
likewise considered in this analysis. 

The leading D — 2 corrections induced by the hght-quark masses are computed using 
the running quark masses evaluated at the two-loop level. The following values are used 
for the renormalization- invariant quark mass parameters TJii defined in Ref. [67]: 

m„ = (8.7 ± 1.5) MeV , ma = (15.4 ± 1.5) MeV , = (270 ± 30) MeV . (72) 

The dimension D — A operators have dynamical contributions from the gluon 
condensate {{as/7i)GG) and quark condensates mu{0\uu\0), m(i(0|(i(i|0) of the light u,d 
quarks. The contribution of the gluon condensate to Rr,v/A vanishes in first order in 
as{ml). However, there appear second order terms in the Wilson coefficients due to 
the logarithmic s dependence of as{s) which after performing the integral (66) becomes 
Q;^(m^). The quark condensates are parametrized [67] by (rrij'ipiipi) — —mjji\ with 

fi^^jl^^ (189 ± 7) MeV , /t^ = (160 ± 10) MeV . (73) 

The contributions from dimension D = & and 8 operators are rather complex. They 
are treated phenomenologically in this analysis. Higher order contributions from D > 10 
operators are expected to be small as, equivalent to the gluon condensate, constant terms 
and terms in leading order vanish in Eq. (66) after integration. 

The formulae are taken entirely from Ref. [67], in which Eq. (66) is evaluated after 
the power terms (68) are inserted into the integral. 



17.2 Spectral moments 

It was shown in Ref. [88] that it is possible to benefit from the information provided by 
the explicit shape of the spectral functions in order to obtain additional constraints on 
as(so) and — more importantly — on the nonperturbative condensates. The spectral 
moments at are defined as: 

-"^'^-h(^-^)(^y-^^ (74) 

with -R^V/A ~ Rt,v/a- The factor (1 — s/m^)^ suppresses the integrand at the crossing of 
the positive real axis where the validity of the OPE is less certain and the experimental 
accuracy is statistically limited. Its counterpart {s/m^Y projects out higher energies. 
The spectral information is used to fit simultaneously as{m'^) and the phenomenological 
effective operators {{ols/t^)GGd=a) i {Od^q) and {Od=^). Due to the intrinsic strong 
correlations only five moments are used as input to the fit. 

In analogy to i?,-, the contributions to the moments originating from perturbative 
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and nonperturbative QCD are separated via the OPE. The prediction of the perturbative 
contribution takes the form 

3 

Sm) = ^ir„A("''=')(Q;,) , (75) 

n=l 

with contour integrals A^'^'''''\as) [88] that are expanded up to a^{s) (FOPT) or 
numerically resolved for the running as{—s) obtained from the RGE (FOPTci). 

For practical purpose it is more convenient to define moments that are normalized to 
the corresponding Rt,v/a in order to decouple the normalization from the shape of the r 
spectral functions: 

^.'u - 1^ - "h-^ i^-^M ^\ . (76) 



' Rr.viA I \ m?J \mij Nv,A ds 

The two sets of experimentally almost uncorrelated observables — Rt,v/a on one hand 
and the spectral moments on the other hand — yield independent constraints on asirri^) 
and thus provide an important test of consistency. 

The results given in Section 14.5 are recalled: 

Rr,v = 1-787 ± 0.013 , (77) 
Rr,A = 1-695 ± 0.013 , (78) 
R^y+A = 3.482 ± 0.014 . (79) 

The normalization according to Eq. (76) reduces considerably the correlations between 
Rr and the moments. It is completely neghgible in the {V + A) case where Rr,v+A is 
calculated from the difference Rt — Rt,s, which has no correlation with the hadronic 
invariant mass spectrum. 

The measured values of the moments for the V, A and {V + A) spectral functions are 
given in Table 23, and their correlation matrices in Table 24. The correlations between 
the moments are computed analytically from the contraction of the derivatives of two 
involved moments with the covariance matrices of the respective normalized invariant 
mass-squared spectra. In all cases, the negative sign for the correlations between the 
k = 1, / = and the A; = 1, / > 1 moments is understood to be due to the p and the vr, 
Oi peaks which determine the major part of the k = 1,1 = moments. They are much 
less important for higher moments and consequently the amount of negative correlation 
increases with I — 1,2,3. This also explains the large and increasing positive correlations 
between the k — 1,1 > 1 moments, in which, with growing I, the high energy tail becomes 
more important than the low energy peaks. The total errors for the {V + A) case are 
dominated by the uncertainties on the hadronic fractions. 



17.3 Results of the fits: ^^(m^) 

The combined fits to the measured V, A and {V + A) ratios Rj- and moments adjust 
the parameters Q;s(m^), {{as/7r)GG), {Oq)v/a and {Os)v/a of the OPE in the theoretical 
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Table 23: Spectral Moments of vector (V), axial- vector (A) and vector plus axial- vector 
{V + A) inclusive r decays. The errors give the total experimental uncertainties including 
statistical and systematic effects. 



ALEPH 


/ = / = ! J = 2 1 = 3 




0.7160 0.1694 0.0533 0.0227 
0.0025 0.0004 0.0006 0.0004 


r)i/ 


0.7107 0.1489 0.0657 0.0317 
0.0018 0.0006 0.0004 0.0003 


^V+A 


0.7134 0.1595 0.0594 0.0271 
0.0016 0.0004 0.0003 0.0003 



Table 24: Experimental correlations between the moments D^y/A/v+A- There are no 
correlations between Rry+A and the corresponding moments. 



ALEPH 


r)10 

^ry 


nil 
^ry 


r)i2 
^ry 


r)i3 
^ry 


ALEPH 


Dl^A 




^5 


r)i3 


R-T-y 


-0.34 


0.20 


0.33 


0.35 


Rt,A 


-0.23 


0.02 


0.16 


0.25 


r)io 
^ry 


1 


-0.82 


-0.98 


-0.99 


r)io 


1 


-0.78 


-0.97 


-0.98 


^ry 




1 


0.90 


0.82 


nil 




1 


0.87 


0.66 


7-)12 

^ry 
^ry 






1 


0.99 
1 


DI'a 






1 


0.94 
1 



ALEPH 


nlO nil r)12 r)13 
^ry+A ^ry+A ^ry+A ^ry+A 


r)io 

nil 

r)i3 


1 -0.73 -0.98 -0.99 
1 0.80 0.67 

- - 1 0.98 

- - - 1 



predictions (67) and (74) of the above quantities. 

The uncertainties entering the theoretical predictions have been estimated in great 
detail in the previous analysis [12]. Minor contributions originate from the CKM matrix 
element \Vud\, the clcctroweak radiative correction factor Sew, the light quark masses 
mu,d, and the quark condensates. The largest contributions have their origin in the 
truncation of the perturbative expansion. Although it might introduce double-counting, 
the procedure used considers separate variations of the unknown higher order coefficient 
K4 and the choice of the renormalization scale. Thus the estimate is expected to be 
conservative. The renormalization scale is varied around rrir from 1.1 to 2.3 GeV with 
the variation over half of the range taken as systematic uncertainty. Concerning the 
coefficient K4, some theoretical progress has been achieved for its direct calculation. Part 
of the contributions have been recently calculated [94], lending support to the estimate 
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(i^4 ^ 27.5) using the principle of minimal sensitivity [95], which allows the computation 
of a renormalization scheme with optimal convergence [96] . Another estimate of K4 was 
obtained [97] using the a priori freedom of the choice of the renormalization scale // to 
increase the sensitivity of the perturbative series on K4, yielding K4 = 27 ± 5. Taking 
advantage of the new theoretical developments, the value K4 = 25 ± 25 is assumed in this 
analysis, instead of 50 ± 50 used in the previous analysis [12]. 

The fit minimizes the of differences between measured and fitted quantities 
contracted with the inverse of the sum of the experimental and theoretical covariance 
matrices. 

The results are listed in Table 25, separately for the FOPT and FOPTci approaches. 
It is to be noted that even if the fits are better with the FOPT expansion, the two sets 
of results are considered on the same footing and averaged later on. Table 26 gives the 
corresponding correlation matrices between the fitted parameters. The limited number 
of observables and the strong correlations between the spectral moments explain the 
large correlations observed, especially between the fitted nonperturbative operators. The 
precision of as(m^) obtained with the two perturbative methods employed is comparable, 
however their central values differ by 0.01-0.02. The S^"^^ term is the pure theoretical 
contribution from the known masses (within errors) of the light u, d quarks. In the 6^^^ 
term the quark condensates and the quark masses are rather well known and are fixed 
theoretically within their errors, while the gluon condensate is adjusted in the fit. 

As mentioned before, there exists no constraining prescription which allows a 
resolution of the ambiguity between FOPTci and FOPT. The final result on as(m^) 
is thus taken as the arithmetic average of the two values given in Table 25, with half of 
their difference added as additional theoretical error: 

a,(m^) = 0.340 ± 0.005exp ± 0.014th . (80) 

The first error accounts for the experimental uncertainty, the second number gives the 
uncertainty of the theoretical prediction of Rr and the spectral moments as well as the 
ambiguity of the theoretical approaches employed. 

In Table 25, a remarkable agreement within statistical errors is found between the 
as(m^) values using vector and axial- vector data, providing a first important consistency 
check of the analysis, since the two corresponding spectral functions are independently 
measured and manifest a quite different resonant behaviour. The main results are 
displayed in Fig. 72. 

17.4 Results of the fits: nonperturbative contributions 

The advantage of separating the vector and axial-vector channels and comparing to the 
inclusive {V + A) fit becomes obvious in the adjustment of the leading nonperturbative 
contributions of dimension D = 6 and D = 8, which approximately cancel in the inclusive 
sum. This cancellation of the nonperturbative terms increases the confidence in the 
as{ml) determination from the inclusive {V + A) observables. The gluon condensate is 
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Table 25: Fit results of Q;s(^r) ^^e OPE nonperturbative contributions from vector, 
axial- vector and {V + A) combined fits using the corresponding ratios Rr and the spectral 
moments as input parameters. Where two errors are given they denote experimental 
(first number) and theoretical uncertainties (second number). The S^'^^ term is the 
pure theoretical prediction with quark masses varying within their prescribed range (see 
text). The quark condensates in the 5^'^^ term are fixed to their theoretical values [12] 
within errors and only the gluon condensate is varied as a free parameter. The total 
nonperturbative contribution is the sum 5np = 5^'^^+5^^^+5^^\ Results are given separately 
for the FOPTci and the FOPT perturbative expansion (see text) . 



fitted variable 


Vector {V) 


Axial- Vector (A) 


V + A 


asimi) (FOPTci) 
as{ml) (FOPT) 


0.355 ± 0.008 ± 0.009 
0.331 ±0.006 ±0.012 


0.333 ± 0.009 ± 0.009 
0.327 ±0.007 ±0.012 


0.350 ± 0.005 ± 0.009 
0.331 ± 0.004 ± 0.012 


i^GG) (GeY') (FOPTci) 
i^GG) (GeV^) (FOPT) 


(0.4 ± 0.3) X 10-^ 
(1.5 ±0.3) X 10-2 


(-1.3 ±0.4) X 10-2 
(-0.2 ± 0.4) X 10-2 


(-0.5 ± 0.3) X 10-2 
(0.6 ± 0.2) X 10-2 


5(6) (FOPTci) 
5(6) (FOPT) 


(2.85 ± 0.22) X 10-2 
(2.70 ±0.25) X 10^2 


(-3.23 ± 0.26) X 10-2 
(-2.96 ±0.31) X 10-2 


(-2.1 ±2.2) X 10-3 
(-1.6 ±2.5) X 10-3 


(FOPTci) 
5(8) (FOPT) 


(-9.0 ±0.5) X 10-3 
(-8.6 ±0.6) X 10-3 


(8.9 ±0.6) X 10-3 
(8.6 ± 0.6) X 10-3 


(-0.3 ±4.8) X lO-'' 
(1.2 ±5.2) X 10-4 


xVl DF (FOPTci) 
xVl DF (FOPT) 


0.52 
0.01 


5.71 
0.63 


3.84 
0.11 


^(2) (FOPTci) 
5(2) (FOPT) 


(-3.3 ±3.0) X 10-^ 
(-3.0 ±3.0) X 10-^ 


(-5.1 ±3.0) X lO-'' 
(-5.0 ±3.0) X 10-"^ 


(-4.4 ±2.0) X 10-4 
(-4.0 ±2.0) X 10-4 


SW (FOPTci) 
5(4) (FOPT) 


(4.1 ± 1.2) X 10-"^ 
(6.8 ± 1.0) X 10-4 


(-5.7 ±0.1) X 10-3 
(-5.3 ±0.1) X 10-3 


(-2.7 ±0.1) X 10-3 
(-2.4 ±0.1) X 10-3 


Total 5np (FOPTci) 
Total 5np (FOPT) 


(1.99 ±0.27) X 10-2 
(1.91 ±0.31) X 10-2 


(-2.91 ±0.20) X 10-2 
(-2.63 ±0.25) X 10-2 


(-4.8 ± 1.7) X 10-3 
(-3.9 ±2.0) X 10-3 



Table 26: Correlation matrices according to the fits presented in Table 25 for vector (left 
table), axial- vector (middle) and (V + A) (right table) using FOPTci- As the gluon 
condensate contributes only insignificantly to 5^^\ the correlations to the total 5^4) term 
are small. 



ALEPH 


{GG)v 5^^^ 5f 


{GG)a 5f 5f 


{GG)v+A ^v+A ^V+A 


asim'i) 

{GG)v/A/V+A 
.(6) 

"V/A/V+A 
AS) 

"V/A/V+A 


-0.39 -0.28 -0.34 
1 0.44 0.46 
1 -0.98 
- - 1 


-0.57 0.52 -0.55 
1 -0.81 0.80 
1 -0.99 
- - 1 


-0.37 0.38 -0.45 
1 -0.65 0.65 
1 -0.98 
1 
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0.3 0.32 0.34 0.36 0.38 0.4 

Figure 72: Results for ^^(m^) from the fits of RTy,A,v+A and the moments -Dy^ using 
the FOPTciand FOPT perturbative expansions. The measurements are correlated due 
to the theoretical errors (see Table 25). 

determined by the first / = 0, 1 moments, which receive lowest order contributions. It 
is observed that the values obtained in the V and A fits are not very consistent, which 
could indicate some problem in the validity of the OPE approach, since nonperturbative 
contributions are significant in these cases. Taking the value obtained in the {V + A) 
inclusive fit, where nonperturbative effects are negligible, and adding as systematic 
uncertainties half of the difference between the vector and axial-vector fits as well as 
between the FOPTci and FOPT results, the gluon condensate is found to be 

{{aj7r)GG) = (0.001 ± 0.012) GeV^ (81) 

This result does not provide evidence for a nonzero gluon condensate, but it is consistent 
with and has comparable accuracy to the independent value obtained using charmonium 
sum rules and e+e^ data in the charm region, (0.011 ± 0.009) GeV^ in a combined 
determination with the c quark mass [102]. 

The D = 6, 8 nonperturbative contributions are obtained after averaging the FOPT 
and FOPTci values: 

S^v^ = (2.8 ± 0.3) X 10-2 , S^v^ = (-8.8 ± 0.6) x 10"=^ , 
5f = (-3.1 ± 0.3) X 10-2 , 5f = (8.7±0.6) X 10-^ (82) 
4+A = (-1.8 ± 2.4) X 10-=^, 4+A = (0.5 ± 5.1) X 10-4 . 

The remarkable feature is the approximate cancellation of these contributions in the V + A 
case, both for D = 6 and D — 8. This property was predicted [67] ior D — 6 using 
the simplifying assumption of vacuum saturation for the matrix elements of four-quark 
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Figure 73: The ratio Rry+A versus the square "r mass" Sq. The curves are plotted as 
error bands to emphasize their strong point-to-point correlations in sq. Also shown is the 
theoretical prediction using FOPTci and the results for R^y+A and the nonperturbative 
terms from Table 25. 

operators, yieding 6y^ = —(7/11) 6^^^ and Sy'^ = (2.5 ±1.3) x 10~^, in fair agreement with 
our results. There is no reliable prediction for the D = 8 contributions. 

The total nonperturbative V + A correction is 5np,v+a = (—4.3 ± 1.9) x 10~^, 
much smaller than the corresponding values in the V and A components, ^NP.y = 
(2.0 ±0.3) X 10-^ and ^Npy = (-2.8 ±0.3) x IQ-'^. After removing the theoretically well- 
known contribution from the quark condensates {D = 4), the remaining non perturbative 
component in ^ ± A is —(1.7 ± 1.9) x 10~^, one order of magnitude smaller compared to 
V and A, and even consistent with zero. 

17.5 Evolution of ag with the mass scale 
17.5.1 Test of the running of as{s) at low energies 

Using the spectral functions, one can simulate the physics of a hypothetical r lepton 
with a mass y/so smaller than rrir through Eq. (66). Assuming quark-hadron duality, 
the evolution of Rt{so) provides a direct test of the running of as{so), governed by the 
RGE /5-function. On the other hand, it is a test of the validity of the OPE approach in 
T decays. The studies performed in this section employ only FOPTci. Results obtained 
with FOPT are similar and differ only in the central as(m^) value. 

The functional dependence of Rry+Aiso) is plotted in Fig. 73 together with the 
theoretical prediction using the results of Table 25. The slight mismatch appearing at 
So = rn^ between the experimental and theoretical values of Rry+A is due to the fact that 
the global fits presented in Table 25 involve not only the measured Rry+A, but also the 
spectral moments. The spread due to uncertainties are shown as bands. The correlations 
between two adjacent bins si < S2 are large as the only new information is provided 
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Figure 74: The running of as(so) obtained from the fit of the theoretical prediction to 
Rt,v+a{so). The shaded band shows the data including experimental errors. The curve 
gives the four-loop RGE evolution for three flavours. 

by the small mass difference between the two bins and the slightly different weight 
function. They are reinforced by the original experimental and theoretical correlations. 
Below 1 GeV^ the error of the theoretical prediction of Rt,v+a{so) starts to blow up 
due to the increasing uncertainty from the unknown perturbative term; errors of the 
nonperturbative contributions are not contained in the theoretical error band. Figure 74 
shows the plot corresponding to Fig. 73, translated into the running of ^^(so), i.e., the 
experimental value for as(so) has been individually determined at every sq from the 
comparison of data and theory. Also plotted is the four-loop RGE evolution using three 
quark flavours. 

It is remarkable that the theoretical prediction using the parameters determined at 
the r mass and Rry+A{so) extracted from the measured V + A spectral function agree 
down to So ~ 0.8 GeV^. The agreement is good to about 2% at 1 GeV^. This result, 
even more directly illustrated by Fig. 74, demonstrates the validity of the perturbative 
approach down to masses around 1 GeV, well below the r mass scale. The agreement with 
the expected scale evolution between 1 and 1.8 GeV is an interesting result, considering 
the relatively low mass range, where is seen to decrease by a factor of 1.6 and reaches 
rather large values ~ 0.55 at the lowest masses. This behaviour provides confidence that 
the asi'fn'^) measurement is on solid phenomenological ground. 

17.5.2 Evolution to Mz 

The evolution of the as(m^) measurement from the inclusive {V + A) observables to the 
Z mass scale can be performed, based on the Runge-Kutta integration of the differential 
equation to four loops and matching conditions at heavy quark (c, b) thresholds to three 
loops [98, 100, 101, 103, 104]. There is some ambiguity in the procedure followed to 
satisfy the matching conditions, in particular whether the matching should be performed 




132 



at a scale given by the quark mass or twice the quark mass, and also on the choice of 
quark mass, i.e. the running mass or the pole mass. Different procedures have been tried 
with the relevant parameters varied within their uncertainties. The results are consistent 
within the quoted systematic uncertainty and for definiteness the final value is obtained 
with the matching realized at the quark-pair thresholds: 

a^(M|) = 0.1209 ± 0.0006exp ± 0.0016th ± O.OOOSevoi 

= 0.1209 ± 0.0018 . (83) 

The first two errors originate from the asim^) determination given in Eq. (80), while 
the last error stands for possible ambiguities in the evolution due to uncertainties in 
the matching scales of the quark thresholds. The uncertainty on the evolution has 
main components from the c-quark mass (0.00020), the 6-quark mass (0.00005), the 
renormalization scale (taken from 0.7 to 3.0 nig, 0.00023), the 3-loop truncation in 
the matching expansion (0.00026) and the four-loop truncation in the RGE equation 
(0.00031). 

The result (83) is a 1.5% determination of the strong couphng constant at the Mz 
scale, limited in accuracy by theoretical uncertainties in the pcrturbative expansion. 
The significant improvement in precision compared to our previous determination [12], 
0.1202 ± 0.008exp ± 0.0024th ± O.OOlOevoi, results from the higher statistics and the more 
detailed experimental analysis, but mostly from the increased knowledge from theory, 
essentially at the level of the perturbative expansion [94] , and also for the scale evolution 
procedure. 

The q;s(M|) result from the present analysis of r decays is in agreement with the 
determination from the Z width as obtained by the global fit to all electroweak data [105]: 

«s(M|)z width = 0.1186 ± 0.0027 , (84) 

as demonstrated in Fig. 75. Both determinations have somewhat comparable accuracies, 
which are however very different in nature. On one hand, the r value is dominated 
by theoretical uncertainties. On the other hand, the determination at the Z resonance, 
benefiting from the much larger energy scale and the correspondingly small uncertainties 
from the truncated perturbative expansion, is limited by the experimental precision on the 
electroweak observables, essentially the ratio of leptonic to hadronic peak cross sections. 
The consistency between the two results, Eqs. (83) and (84), provides a powerful test, over 
a range of s spanning more than three orders of magnitude, of the evolution of the strong 
interaction coupling as predicted by the nonabelian nature of the QCD gauge theory. It 
is therefore a precise test of asymptotic freedom at the 2.7% level. 

18 Conclusions 

A final analysis of r decay branching fractions has been presented using all LEP 1 data 
collected with the ALEPH detector. As in the pubhcation based on the 1991-1993 data set. 
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Figure 75: The evolution of the measured as(m^) to higher scales fi using the 4-loop RGE 
and the 3-loop matching conditions applied at the heavy quark-pair thresholds (hence the 
discontinuities at 2 rric and 2 rrih). The evolution of the ALEPH Q;s(m^) determination is in 
agreement with the direct value obtained from the Z width through the global electroweak 
fit [105]. 

it uses a global analysis of all modes, classified according to charged particle identification, 
and charged particle and 7r° multiplicity up to 4 ir^s in the final state. Major improvements 
are introduced with respect to the published analysis and a better understanding is 
achieved, in particular in the separation between genuine and fake photons. In this process 
shortcomings and small biases of the previous method were discovered and corrected, 
leading to more robust results. As modes with kaons [K^, Kg, and K^) have already 
been studied and pubhshed with the full statistics, the non-strange modes without kaons 
are emphasized in this paper. Taken together these results provide a complete description 
of r decays up to 6 hadrons in the final state. All decay channels with branching fractions 
larger than 0.1% are considered in this analysis. 

The measured branching fractions are internally consistent and agree with known 
constraints from other measurements in the framework of the Standard Model. More 
specifically, the results on the leptonic and pionic fractions lead to powerful tests of 
universality in the charged leptonic weak current, showing that the e — n — t couplings 
are equal within 0.3%. The branching ratio of r — > Ut-tttt^, which is of particular interest to 
the accurate determination of vacuum polarization effects, is determined with a precision 
of 0.5% to be (25.47 ± 0.13) %. Also the ratio of oi branching fractions into 7r27r° and 
37r final states is measured to be 0.979 ± 0.018, in agreement with expectation from 
partial wave analyses of these decays. Separating nonstrange hadronic channels into 
vector (V) and axial-vector (A) components and normalizing to the electronic width 
yields the ratios Rr,v = 1-787 ± 0.013, Rr,A = 1-695 ± 0.013, Rry+A = 3.482 ± 0.014 and 
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R^y^A = 0.092 ± 0.023. 

Spectral functions for the hadronic modes are determined through the unfolding of 
the measured invariant mass distributions, taking into account detector effects. As for 
the branching fractions, the spectral functions are separated into vector and axial-vector 
contributions. The vector spectral functions are compared to the corresponding ones 
obtained through e+e~ annihilation, providing accurate tests of isospin invariance of the 
electroweak current. Both spectral functions are fitted within QCD to pcrturbative and 
nonperturbative components. The latter are determined to be very small, while the 
dominant perturbative part yields a precise value for the strong coupling at the r mass 
scale, as{ml) — 0.340 ± (0.005)ea;p ± (0.014)(ft. The evolution to the Z mass scale yields 
as{M'^) = 0.1209 ± 0.0018. It agrees well with the direct determination from the Z width 
and provides the most accurate test to date of asymptotic freedom in the QCD gauge 
theory. 
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